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ABSTRACT 

Theoretical light curves of four recurrent novae in outburst are modeled to obtain various physical 
parameters. The four objects studied here are those with a red giant companion, i.e., T Coronae Borealis, 
RS Ophiuchi, V745 Scorpii, and V3890 Sagittarii. Our model consists of a very massive white dwarf 
(WD) with an accretion disk and a red giant companion. Light curve calculation includes reflection 
effects of the companion star and the accretion disk together with a shadowing effect on the companion 
by the accretion disk. We also include a radiation-induced warping instability of the accretion disk to 
reproduce the second peak of T CrB outbursts. The early visual light curves are well reproduced with 
a thermonuclear runaway model on a very massive white dwarf close to the Chandrasekhar mass limit, 
i.e., MwD = 1-37 ± 0.01 Af© for T CrB, 1.35 ± 0.01 Mq for RS Oph with solar metallicity {Z = 0.02) 
but 1.377 ± 0.01 Mq for RS Oph with low metallicity {Z = 0.004), 1.35 ± 0.01 Mq for V745 Sco, and 
1.35 ± 0.01 Mq for V3890 Sgr. Optically thick winds, which blow from the WDs during the outbursts, 
play a key role in determining the nova duration and the speed of decline because the wind quickly 
reduces the envelope mass on the WD. Each envelope mass at the optical maximum is also estimated 
to be AM - 3 X IO-'^Mq (T CrB), 2 x IO-'^Mq (RS Oph), 5 x IO-^Mq (V745 Sco), 3 x 10'^ Mq 
(V3890 Sgr), indicating an average mass accretion rate Mace ^ 0.4 x 10~'^Mq yr~^ (80 yrs, T CrB), 
1.2 X IQ-'^Mq yr-i (18 yrs, RS Oph), 0.9 x 10-^Mq yr"! (52 yrs, Y745 Sco), 1.1 x IO'^Mq yr"! (28 yrs, 
V3890 Sgr) during the quiescent phase before the last outburst. Although a large part of the envelope 
mass is blown in the wind, each WD can retain a substantial part of the envelope mass after hydrogen 
burning ends. Thus, we have obtained net mass-increasing rates of the WDs as A/rc 0.1 x lO-'^M© 
yr-i (T CrB), 0.12 x lO-'^Mg yr'^ (RS Oph), 0.05 x lO-^M© yr"! (V745 Sco), 0.11 x IQ-'^Mq yr"! 
(V3890 Sgr). These results strongly indicate that the WDs in the recurrent novae have now grown up 
to near the Chandrasekhar mass limit and will soon explode as a Type la supernova if the WDs consist 
of carbon and oxygen. We have also clarified the reason why only T CrB shows a secondary maximum 
but the other three systems do not. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: supernovae — stars: 
symbiotic — stars: individual (RS Oph, T CrB, V3890 Sgr, V745 Sco) 



1. INTRODUCTION 

Type la supernovae (SNe la) are one of the most lumi- 
nous explosive events of stars. Recently, SNe la have been 
used as good distance indicators that provide a promising 
tool for determining cosmological parameters (Perlmutter 
et al. 1999; Riess et al. 1998) because of their almost uni- 
form maximum luminosities. These both groups derive the 
maximum luminosities, imax, of SNe la completely empir- 
ically from the light curve shape (LCS) of nearby SNe la, 
and assume that the same Lmax-LCS relation holds in high 
red-shift (z) SNe la. In order that this method works, the 
nature of SNe la and their progenitors should be the same 
between high-z and low-z SNe la. Therefore, it is neces- 
sary to identify the progenitor binary systems, thus con- 
firming whether or not SNe la "evolve" from high red-shift 
to low red-shift galaxies (e.g., Umeda et al. 1999a). 

It is widely accepted that the exploding star itself is 



an accreting white dwarf (WD) in a binary (e.g., Nomoto 
1982; Nomoto, Thielemann, & Yokoi 1984). However, the 
companion star (and thus the observed binary system) is 
not known. Several types of objects have ever been con- 
sidered, which include merging double white dwarfs (e.g., 
Iben & Tutukov 1984; Webbink 1984), recurrent novae 
(e.g., Starrfield, Sparks, & Truran 1985), symbiotic stars 
(e.g., Munari & Renzini 1992) etc (see, e.g., Livio 2000, 
for recent summary). 

One of the candidates, the recurrent nova U Scorpii 
underwent sixth recorded outburst in February 25, 1999 
(Schmeer 1999). For the first time, a complete light curve 
of the outburst has been obtained from the rising phase 
to the final decline phase through the mid-plateau phase 
including eclipses (e.g., Hachisu et al. 2000a; Matsumoto, 
Kato, & Hachisu 2000; Munari et al. 1999). By analyz- 
ing these new data on a basis of the theoretical light-curve 
model of outbursts, Hachisu et al. (2000a) have elucidated 
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the outburst nature of U See much more accurately than 
the previous estimates (e.g., Starrfield, Sparks, & Shaviv 
1988; Kato 1990). We summarize some critically impor- 
tant values they have obtained. (1) The white dwarf mass 
in U Sco is l.STitO.OlM©. (2) The companion is a slightly 
evolved main-sequence star of 1.4 — 1.6 Mq that expands 
to fill the Roche lobe after a large part of the central hy- 
drogen has been consumed and is now transferring mass 
to the white dwarf in a thermally unstable manner. (3) 
The white dwarf is now growing in mass at an average rate 
of ~ 1 X 1O~^M0 yr~^. Therefore, the white dwarf will 
reach the critical mass (Mja = I.STSMq, in W7 model of 
Nomoto et al. 1984) for the SN la explosion in the quite 
near future. They concluded that U Sco is a very strong 
candidate for the immediate progenitor of SNe la. Other 
recurrent novae are also expected to be a strong candidate 
for SNe la. 

Thus, we have started our light curve analysis of recur- 
rent novae in which the mass and the growth rate of white 
dwarfs are determined as accurately as possible and ex- 
amine the possibility of an SN la explosion in the quite 
near future. In this series of papers, details of our results 
on recurrent novae are presented, although a part of our 
results has already been published in several short papers 
(Hachisu & Kato 1999 for T CrB; Hachisu et al. 2000a,b 
for U Sco; Hachisu & Kato 2000a for RS Oph; Hachisu & 
Kato 2000b for V394 CrA). 

Here;, in the; first paper of this series, we have examined 
four recurrent novae which have a red giant companion, 
T CrB, RS Oph, V745 Sco, V3890 Sgr. Two of them (T 
CrB and RS Oph) are reexamined with new observational 
knowledge. We first summarize the nature of recurrent no- 
vae and explain the physical reason why they become so 
massive as near the Chandrasekhar mass limit in §2. Then, 
we show details of our numerical methods for calculating 
evolution of bloated white dwarf envelopes in §3 and for 
obtaining theoretical light curves in §4, because details of 
our numerical method have not been fully described yet. 
New numerical results on determination of white dwarf 
masses for four recurrent novae are shown and discussed 
in §5 for T CrB, in §6 for RS Oph, in §7 for V745 Sco, and 
in §8 for V3890 Sgr. Conclusions follow in §9. 

2. RECURRENT NOVAE AS A PROGENITOR OF TYPE lA 

SUPERNOVAE 

Recurrent novae are characterized by nova eruptions 
with its recurrence time scale from a decade to a century, 
its brightness amplitude larger than 6 magnitudes, and its 
expansion velocity faster than 300 km s~^ (e.g., Webbink 
et al. 1987). They are considered to be binary systems in 
which the secondary companion star overfiows the inner 
critical Roche lobe and transfers matter onto the primary 
white dwarf (WD) component. When the hydrogen-rich 
envelope on the WD reaches a critical mass, hydrogen ig- 
nites to trigger a nova eruption (hydrogen shell-flash). We 
call this a thermonuclear runaway (TNR) event. 

Recurrent novae are morphologically divided into three 
groups depending on its companion star: dwarf com- 
panions, only one member, T Pyxidis with P = 0.076 
days, (Schaefer 1990, 1992; Patterson et al. 1998); 
slightly evolved main-sequence (MS) companions, U Scor- 
pii with P = 1.23 days (Schaefer 1990; Schaefer & Ring- 



wald 1995), V394 Coronae Austrinae with P = 0.757 
days (Schaefer 1990), Large Magellanic Cloud Nova 1990 
No.2, hereafter LMC RN, with P =unknown but P ~ 1.3 
days was suggested (Sekiguchi 1992), CI Aquilae with 
P = 0.618 days (Mennickent & Honeycutt 1995), and 
red giant (RG) companions, T Coronae Borealis with 
P = 227.6 days (Lines, Lines, & McFaul 1988), RS Ophi- 
uchi with P = 460 days, (Dobrzycka & Kenyon 1994), 
recently revised to be P = 457 days (Fekel et al. 2000), 
V745 Scorpii with P =unknown, V3890 Sagittarii with 
P =unknown, where P is the orbital period of the binary 
in units of day. These are sometimes called T Pyx sub- 
class, U Sco subclass, and T CrB (or RS Oph) subclass, 
respectively (e.g., Warner 1995). We have summarized 
these features in Table 1. 

Classical novae, sharing many aspects with the recur- 
rent novae, are also binary systems consisting of a white 
dwarf and a less massive main-sequence star that fills the 
Roche lobe and transfers its mass to the white dwarf. The 
white dwarf accumulates hydrogen-rich matter and ignites 
hydrogen atop the white dwarf when the hydrogen-rich en- 
velope on the white dwarf reaches a critical mass. A part 
of the hydrogen diffuses into the white dwarf before the 
ignition so that a very surface layer of the white dwarf 
is dredged up into the hydrogen-rich envelope and blown 
off in the outburst wind (e.g., Prialnik 1986; Kovetz & 
Prialnik 1994). As a result, nova ejecta contain white 
dwarf matter, i.e., c:arbon and oxygen for carbon-oxygen 
white dwarfs or oxygen, neon, and magnesium for 0-Ne- 
Mg white dwarfs (see, e.g., Prialnik 1986; Kovetz &: Pri- 
alnik 1994; Warner 1995, and references therein). Thus, 
the white dwarf is eroded gradually after many cycles of 
nova outbursts, recurrence periods of which are usually 
longer than thousands of years. 

Recurrent novae show some characteristic differences: 
(1) heavy elements such as carbon, oxygen and neon are 
not enriched in ejecta but similar to the solar values (e.g., 
Livio & Truran 1992; Warner 1995, for summary), indi- 
cating that the white dwarf is not eroded; (2) very short 
recurrence periods from a decade to a century theoreti- 
cally require very massive white dwarfs close to the Chan- 
drasekhar mass limit (e.g.. Fig. 9 of Nomoto 1982). These 
two arguments indicate that the mass of the white dwarfs 
in recurrent novae increases toward the Chandrasekhar 
mass limit (e.g., Starrfield, Sparks, & Truran 1985; Star- 
rfield, Sparks, & Shaviv 1988; Livio & Truran 1992: Delia 
Valle & Livio 1996). We expect that, if the white dwarf is 
made of carbon and oxygen, the white dwarf ignites carbon 
at the center and explodes as a Type la supernova when 
it reaches the critical mass of Mja = 1.378-M0 (Nomoto 
1982; Nomoto et al. 1984). However, the evolutionary 
path to recurrent novae has not long been elucidated as 
well as the reason why so massive white dwarfs can exist 
in recurrent novae. 

2.1. Accretion Wind Evolution vs. Common Envelope 

Evolution 

In this subsection, we briefly explain how white dwarfs 
in recurrent novae have grown its mass toward the Chan- 
drasekhar mass limit. The reason why the standard binary 
evolution theory (will be cited below) has failed to explain 
the evolutionary path to recurrent novae is mainly in the 
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theoretical prediction of a common envelope formation and 
the ensuing spiral-in process, because these processes in- 
hibit the growth of white dwarfs in binary systems. As 
discussed in many previous papers on binary evolutions, 
it had been widely accepted that hydrogen-rich envelopes 
on mass-accreting white dwarfs expand to a red giant size 
when its mass accretion rate, Maco exceeds a critical limit, 
Mcr, i.e., A4cc > Afcr ^ 1 X IO-6M0 yr-1 (sec, e.g.. Fig. 
9 of Nomoto 1982), and easily forms a common envelope 
(e.g., Nomoto, Nariai, & Sugimoto 1979). Once a common 
envelope is formed, two stars begin to spiral-in each other 
due to viscous drag; thus, producing a double degenerate 
system (e.g., Iben & Tutukov 1984; Wcbbink 1984). 

However, Hachisu, Kato, & Nomoto (1996) found that 
white dwarfs begin to blow strong winds (t^wind ^ 1000 km 
s""'^ and Mwind ^ 1 x 1O~^M0 yr"-'^) when the mass ac- 
cretion rate exceeds the critical rate, i.e.. Mace > as 
illustrated in Figure 1; thus, preventing the binary from 
collapsing. The binary does not spiral-in but keeps its 
separation. Hydrogen burns steadily and therefore the 
helium layer of the white dwarf can grow at a rate of 
Mhb ~ Mcr- The other transferred matter is blown in 
the wind (A/vvind ~ A^acc — Mcr), where Mwind is the wind 
mass loss rate. We call this accretion wind because it be- 
gins to blow when the accretion rate exceeds the critical 
rate. 

Therefore, as shown in Figure 2, we have to revise Figure 
9 of Nomoto (1982) , which has been widely used in binary 
evolution scenarios as basic processes of mass-accreting 
white dwarfs (e.g., Kahabka & van den Heuvel 1997). The 
most important difference of Figure 2 from old Nomoto's 
(1982) is that white dwarfs grow its mass in a much wider 
parameter region. In our new diagram, the status of a 
mass-accreting white dwarf is specified by the following 
three phases; (1) accretion wind phase (Afacc > Mcr); (2) 
steady shell-burning phase (Mgtd < A/acc < Mcr); and (3) 
intermittent shell-flash phase (Mace < Mgtd), where Mstd 
means the lowest limit of mass accretion rate for steady 
hydrogen shell burning. The new growing region of white 
dwarfs. Mace > Mgtd; is much wider than old Nomoto's 
(1982) narrow strip, M^td < Mace < Mcr- 

2.2. Evolutionary Paths to Recurrent Novae and Type la 

Supernovae 

Based on the new mechanism of Figures 1 and 2, 
Hachisu et al. (1999a, b) have found two paths to SNe la 
(see also, Li &; van den Heuvel 1997). In these two paths, 
white dwarfs can accrete mass continuously from the com- 
panion and grow at a rate of Mcr 1 x 1O~^M0 yr~^. 
Therefore, candidates of SN la progenitors arc systems at 
the final stages on these paths; i.e., one is a supersoft X-ray 
source (SSS) consisting of a white dwarf and a lobe-filling 
more massive main-sequence star (WD-I-MS system, e.g., 
Li & van den Heuvel 1997) and the other is a symbiotic 
star consisting of a white dwarf and a mass-losing red gi- 
ant (WD-I-RG system, e.g., Hachisu et al. 1996). Both 
the systems contribute to main parts of the SN la birth 
rate (Hachisu et al. 1999a, b) (sec also Livio 2000, for a 
recent review). 

Hachisu et al. (1999a, b) have followed many evolution- 
ary paths and obtained the initial parameter regions of 



binaries that finally produce an SN la. Figure 3 shows 
such two regions of the WD-I-MS systems (SSS channel) 
and the WD-I-RG systems (symbiotic channel). We have 
added in this figure the final parameter regions just before 
an SN la explosion occurs. Starting from the initial region 
encircled by the thin solid lines, binary systems evolve and 
can explode as an SN la in the regions enclosed by the thick 
solid lines. Hachisu et al. (1999a,b) showed that some re- 
current novae lie on the boarder of their parameter region 
that can produce an SN la. Among six recurrent novae 
with known orbital periods, five fall in the final regions of 
SN la progenitors just before the explosion. Only T Pyx 
is far out of the regions of SNe la. 

3. OPTICALLY THICK WINDS IN DECAY PHASE OF NOVAE 

In the thermonuclear runaway (TNR) model of nova 
explosions, WD envelopes expand greatly as large as ^ 
100 Rq or more and then the photospheric radius gradu- 
ally shrinks to the original size of the white dwarfs (e.g., 
i?WD = 0.0032 Rq for Mwd = 1-37 Mq). Here, we call 
this decay phase, i.e., the phase after the maximum ex- 
pansion of the photosphere. In the decay phase of novae, 
the photospheric temperature increases from Tph ~ 10^ K 
to ~ 10^ K with the bolomctric luminosity being almost 
constant (near the Eddington limit). Therefore, the main 
emitting region moves from optical to soft X-ray through 
UV. Correspondingly, the optical luminosity reaches its 
maximum at the maximum expansion of the photosphere 
and then decays toward the level in quiescent phase. Since 
the WD envelope reaches a steady-state in the decay phase 
of novae (e.g., Prialnik 1986; Kato & Hachisu 1994), we 
arc able to approximate the development of the envelope 
structure with a unique sequence of steady-state wind so- 
lutions having a different envelope mass (AM) as shown 
by Kato & Hachisu (1994). 

Optically thick winds, blowing from WDs in the decay 
phase of novae, play a key role in determining the nova 
duration, because the winds eject a large part of the enve- 
lope mass and drastically reduce the fuel. The winds blow 
when the WD photosphere is larger than ^ 0.05 — O.li?0, 
i.e., when the photospheric temperature is lower than 
logTph ^ 5.5, because the wind is driven by a strong peak 
at logT 5.2 in OPAL opacity (e.g., Iglesias & Rogers 
1996), and logTph ~ 5.5 corresponds to the shoulder of 
the peak on the high temperature side. This opacity peak 
strongly blocks photons and, as a reaction, matter is ef- 
fectively accelerated to blow a wind deep inside the pho- 
tosphere (see, e.g., Figs. 2, 3, and 4 in Kato & Hachisu 
1994). 

Assuming the solar abundance of heavy elements (Z = 
0.02), we have calculated the steady-state envelope solu- 
tions and obtained sequences for WDs with various masses 
of Mwd = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.34, 1.35, 
1.36, 1.37, and 1.377M0. Here, we have adopted the up- 
dated OPAL opacity (Iglesias & Rogers 1996) with lower 
hydrogen contents oi X = 0.70, 0.50, 0.35, 0.15, 0.10, 0.08, 
0.07, 0.06, 0.05, and 0.04. The numerical methods and var- 
ious assumptions are essentially the same as those in Kato 
& Hachisu (1994) except the opacities. Some of them are 
plotted in Figure 4. 

For RS Oph and LMC RN, we have also obtained the 
envelope solutions for lower metallicities, i.e., Z = 0.001, 
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0.002, 0.004, and 0.01 with X = 0.70, 0.50, 0.35, 0.15, 
0.10, 0.08, 0.07, 0.06, 0.05, and 0.04. Some of them are 
also shown in Figure 5. 

3.1. Evolution of White Dwarf Envelope 

Hydrogen shell-flashes occur when the mass transfer rate 
is below the lowest limit for the steady hydrogen shell- 
burning, i.e.. Mace < Msid- In this subsection, we describe 
such a rapid development of white dwarf envelopes. As 
already mentioned above, optically thick winds blow in 
the decay phase of hydrogen shell-flashes. Since each wind 
solution is a unique function of the envelope mass AM 
for a given WD mass, we can follow the development of 
the envelope structure by calculating the envelope mass, 
which is decreasing due to the wind mass loss at a rate 
of Mwind(AM) and hydrogen shell burning at a rate of 
Mnuc(AM), i.e., 

^AM = Mace - Mwind - Mnue, (1) 

where Mace is the mass accretion rate of the WD. Directly 

integrating equation (1), we obtain the time-development 
of the envelope, AM, and then the evolutionary changes of 
physical quantities such as the photospheric temperature 
Tph(AM), photospheric radius i?ph(AM), photospheric 
velocity t;ph(AM), wind mass loss rate Mwind(AM), and 
nuclear burning rate Mnuc(AM), all of which are a unique 
function of the envelope mass, AM. 

We have shown, in Figures 4 and 5, the envelope mass 
AM, the photospheric radius i?ph, the photospheric tem- 
perature Tph, and the photospheric velocity Vph against the 

mass decreasing rate of the envelope, Mwind + Mnue- The 
wind mass loss rate becomes as large as 10^'* 1O~^M0 
yr~^ at the maximum expansion of the photosphere i?ph ~ 
lOOi?0, and then decreases down to 1 x 1O~^M0 yr~^ just 
before the wind stops, i.e., when i?ph ~ O.li?0. The pho- 
tospheric temperature increases from 1 x 10^ K to 3 x 10^ 
K. The photospheric velocity remains as high as ~ 1000 
km s~^ during the wind phase. 

When the envelope mass decreases to below the critical 
mass, the wind stops and the photosphere rapidly shrinks 
towards its original size before the outburst, for exam- 
ple, from _Rph = O.O6-R0 to O.OO5i?0 within a week for 
-^WD = 1-37 Mq. Shortly after the wind stops, steady- 
state shell burning on the WD ends but hydrogen burning 
itself is still continuing and supplying a part of the lu- 
minosity. In this phase, the envelope mass is decreased 
only by nuclear burning. The rest of the luminosity comes 
from the thermal energy of the hydrogen envelope and the 
hot ash (helium) below the hydrogen layer. The thermal 
energy amounts to several times 10'*"' ergs for the case of 
MwD = 1.37 M0, X = 0.05, and Z = 0.02 (for the 1999 
outburst of U Sco). It can supply a bolometric luminosity 
of 10^^ ergs s~^ for about ten days or so. 

The mass lost by the wind, AM^ind, and the mass added 
to the helium layer of the WD, AMne, are calculated from 

AMwind = j Mwind dt, (2) 

and 

AMHe = / Mnue dt, (3) 



respectively. The efficiency of hydrogen shell-flashes (nova 
outbursts) is defined by 

^ AMHe _ AMne 

AM^ax ~ AMwind + AMHe ' ^ ' 

where AMmax is the envelope mass at the optical maxi- 
mum. 

Assuming a high constant rate of Afacc ^ 10^'^ Af© yr^^ 
in equation (1), we start the calculation. The envelope 
mass is quickly increased and the optical luminosity soon 
reaches the observational maximum in a day or so. Then, 
we drop the rate to Mace ~ 1 x 1O~^M0 yr~^ or zero and 
further follow the decay phase of nova. Therefore, our ris- 
ing phase is not a realistic one but just a rough sketch. 
Once the envelope has reached the optical maximum, our 
theoretical light curve in the decay phase represents a real- 
istic one because the flow in the envelope is approximately 
in a steady-state as already mentioned above. 

3.2. White Dwarf Envelope in Steady Mass Accretion 

If the mass accretion rate (Mace) does not vary, the enve- 
lope solution is obtained from d{AM)/dt = in equation 
(1). The mass accretion rate is in a balance with the mass 
decreasing rate of the envelope. 

Mace = Mwind + Mnue- (5) 

If the steady mass-accretion rate is given, we can directly 
(without integration) obtain the corresponding state of the 
envelope from Figures 4 and 5. The critical envelope-mass 
decreasing rate, i.e., at the cease point of wind, exactly co- 
incides with the critical mass-accretion rate, i.e., at the be- 
ginning point of wind, Mer = (Mace)er = (Mnue)er, which 

can be approximately expressed by 

Mer = 5.3 X 10-'^-^ - 0.40) Me y,-\ (6) 

for various WD masses (Mwd) and hydrogen contents {X) 
with solar metallicity {Z = 0.02). This relation is reduced 
from our wind solutions for X = 0.70, X = 0.50, X = 0.35, 
X = 0.15, and X = 0.10 {Z = 0.02) and has accuracy 
within to 3 — 5%. 

In the same steady mass-accretion mentioned above, we 
also obtain the lower limit of the mass accretion rate for 
steady hydrogen shell-burning as 

M.td = 2.6 X W-'^-^ - 0.40) M0 yr'^ 

» 0.5Mcr, (7) 

for solar metallicity {Z = 0.02). This relation is reduced 
from our steady-state solutions for X = 0.70, X = 0.50, 
X = 0.35, X = 0.15, and X = 0.10 {Z = 0.02) and has 
accuracy better than 10%. 

4. THEORETICAL LIGHT CURVES 

One of our binary models is graphically shown in Figure 
6. A circular orbit is assumed. Our light curves include the 
contributions of three components: photosphere of WD, 
red giant (RG) companion, and accretion disk (ACDK or 
DK). 
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4.1. White Dwarf Photosphere 

We assume a black-body photosphere of the WD enve- 
lope and calculate the ^-magnitude with a window func- 
tion given by Allen (1973). The photospheric surface is 
divided into 32 pieces in the latitudinal angle (A^ = 7r/32) 
and into 64 pieces in the longitudinal angle {A(j) = 27r/64) 
as illustrated in Figure 6. The contributions to the V light 
are summed up from each piece. We do not consider the 
limb-darkening effect because it is as small as AV < 0.2 
mag. 

We have assumed the luminosity of the WD given by 



1 GMwpMa, 
ivwD = i^wD,o + 7; 5 



(8) 



where the first term is the intrinsic luminosity of the WD 

calculated by our wind and steady hydrogen shell burning 
solutions and the second term is the accretion luminosity 
(e.g., Starrfield, Sparks, & Shaviv 1988); Rwd is the ra- 
dius of the WD (e.g., i?wD = 0.0032i?e for 1.37M0 WD). 
The accretion luminosity is as large as 700 Lq for a mass 
accretion rate of Mace ~ 1 x lO^^M© yr"-'^ onto a 1.37M0 
WD. Therefore, the accretion luminosity itself is negligi- 
bly small during the steady hydrogen shell burning phase 
(iwD.o ^ lO^i©) but may be important after the hydro- 
gen burning diminishes. In such a case, we recalculate the 
surface temperature of the white dwarf, i.e., 

• (9) 



including the accretion luminosity. 

4.2. Companion's Irradiated Photosphere 

We have to include the contribution of the companion 
star when it is strongly irradiated by the WD photosphere. 

The companion star is assumed to fill the inner critical 
Roche lobe as shown in Figure 6 for T CrB, but underfill 
for RS Oph, V745 Sco, and V3890 Sgr as described later. 
Dividing the latitudinal angle into 32 pieces {A9 = 7r/32) 
and the longitudinal angle into 64 pieces ( A(t> = 2n/64), we 
have also summed up the contribution from each area, but 
we neglect both the limb-darkening effect and the gravity- 
darkening effect of the companion star because the both ef- 
fects are negligibly small (AV < 0.2 mag) for the strongly 
irradiated companion. Here, we assume that 50% of the 
absorbed energy is reemitted from the hemisphere of the 
companion with a black-body spectrum at a local temper- 
ature, that is, the efficiency of irradiation is t/rg = 0.5 
unless otherwise specified. 

The original (non-irradiated) photospheric temperature 
of the companion star is obtained from the light curve 
fitting, for example, TRG.org = 3100 K for T CrB. The 
irradiated surface temperature is roughly estimated by 

Vug ^^^2 S + '^T^c,ors, (10) 

where r is the distance from the WD and cos is the incli- 
nation angle of the surface. The irradiated temperature is 

as high as Trg ^ 5000 K for cos d = 0.5, Mwd = 1-37 Mq 
(LwD = 2 X 1038 ergs s'^), and r - 150i?o (T CrB). 
Therefore, the effect of the irradiation becomes very im- 
portant in the late phase of outbursts when the irradiated 
hemisphere faces towards the Earth. 



If the accretion disk around the WD blocks the light 
from the WD photosphere, it makes a shadow on the sur- 
face of the companion star. Such an effect is also included 
in our calculation as explained below. 

4.3. Accretion Disk Surface 

We have to include the luminosity coming from the ac- 
cretion disk when it is irradiated by the WD photosphere. 
Here, we assume a circular accretion disk, the edge of 
which is defined by 

i?disk = aii^ (11) 

where a is a numerical factor indicating a size of the accre- 
tion disk, and Rl the effective radius of the inner critical 
Roche lobe given by Eggleton's (1983) formula. 

The surface of the accretion disk is strongly irradiated 
and dragged outward by the strong wind from the WD 
during the nova outburst. Its thickness may be increasing 
outward by both the ablation and the Kelvin-Helmholtz 
instability. We here approximate such a surface by 



h = /3-Rdisk 



-Rdisk 



(12) 



where h is the height of the surface from the equatorial 
plane, zu the distance on the equatorial plane from the 
center of the WD, and /3 is a numerical factor showing the 
degree of thickness. For comparison, we have calculated 
the three cases of jy, i.e., ly = 9/8 for the case of stan- 
dard accretion disk model by Shakura & Sunyaev (1973), 
ly = 2 and v = 3 ior the case of flaring-up disk (Schandl, 
Meyer-Hofmeister, & Meyer 1997). 

The surface of the accretion disk is divided into 32 pieces 
logarithmically evenly in the radial direction and into 64 
pieces evenly in the azimuthal angle as shown in Figure 6. 
The outer edge of the accretion disk is also divided into 64 
pieces in the azimuthal direction and 2 (or 8) pieces in the 
vertical direction by rectangles. When the photosphere 
of the WD becomes very small, e.g., i?disk/^^ph > 10, we 
attribute the first 16 meshes to the outer region (from 

vj = iidisk to zu = i?disk/'\/lO) to avoid coarse meshes in 
the outer part, and then 16 meshes to the inner region 
(from w = i?disk/"\/TO to tz7 = Rph), each region of which 
is divided logarithmically evenly. 

The surface of the accretion disk absorbs photons and 
reemits in the same way as the companion does. Each area 
of the heated disk surface emits blackbody photons at its 
local temperature, Tdisk- We assume that 50% of the ab- 
sorbed energy is emitted from the surface while the other 
is carried into interior of the accretion disk and eventually 
brought into the WD, that is, the efficiency of irradiation 
is rjuK = 0.5. 

The non-irradiated temperature of the disk surface is 
assumed to be determined by the viscous heating of the 

standard accretion disk model (Shakura & Sunyaev 1973). 
Then, the original disk surface temperature is given by 



<jt: 



SGMwdM^ 



disk, org 



(13) 



The viscous temperature is estimated to be Tdisk, org ~ 
9,000 K at = IRq for Mwd = 1.37M0 and Mace = 
1 X lO~^A'/0 yr^^. However, the viscous heating is neg- 
ligibly small when the disk surface is strongly irradiated 
because of Tdisk ~ 40,000—50,000 K at = IRq for the 
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luminosity of the WD, Lwd = 50, OOOLq, and the average 
inchnatiou of the surface, coaO = 0.1, which is estimated 
from (e.g., Schandl et al. 1997) 



disk 



L\VD f, , SGMwD-^acc 
^Y^V.- ^COS^+ ■ 



, (14) 

The outer rim of the accretion disk is not irradiated by the 
WD photosphere so that the temperature of the disk rim 
is simply assumed, for example, to be 2000 K for T CrB. 

4.4. Radiation-induced Warping 

Radiation-induced warping of accretion disks have been 
suggested by Pringle (1996). Here, we have introduced 
a warping surface calculated by Maloney, Begelman, & 
Pringle (1996). In their formulation, when the tilt vector 
of the accretion disk surface is defined by 

l{wA) = (sin/5cos7, sin^sin7, cos^), (15) 

one of their stationary solutions, 

W = Pe'"', (16) 

can be analytically solved as 

. sin (a;) 



M(l,2,2ix) = e'' 
where x cx So, we have 

and 



f3 = C2 



7 = cittV^, 



(17) 



(18) 



(19) 

(20) 



where ci and C2 are numerical factors indicating a degree 
of warping. Here, we adopt ci = 0.5 and C2 = —1.0, un- 
less otherwise specified. Then, we have transformed the 
coordinates of the disk, i.e., equation (12), by a rotation 
matrix of 

' cos /3 cos 7 — sin 7 — sin (3 cos 7 \ 
cos /3 sin 7 cos 7 — sin /3 sin 7 I . (21) 
sin^ cos ,9 / 

Such an example is shown in Figure 6 for T CrB. As will 
be discussed later, the condition of warping instability is 
satisfied only in T CrB. 

4.5. Numerical Method 

Total V light of the system is calculated from the contri- 
bution of each patch with a window (response) function of 
the F-filter (e.g., Allen 1973). Each patch is almost quad- 
rangular but, strictly speaking, not quadrangular for the 
Roche geometry of the companion star and for the warping 
surface of the accretion disk, because the four vertices are 
not on a plane. Therefore, we define approximately the 
normal unit vector to each patch as follows: the center of 
gravity of each patch can be defined by 



re 



lirA 



(22) 



where G denotes the center of gravity of the patch, and 
A, B, C, and D are the apexes of each quadrangle as illus- 
trated in Figure 7. Defining m as the normal unit vector 
to AABG, n2 to ABCG, 713 to ACDG, and m to ADAG, 



we obtain an approximate normal unit vector of n to patch 
□ABCD by 

_ Smi + S2n2 + SsTis + S4n4 , . 

S,+S2 + S^ + S4 ' ^"^^ 
where each area of the triangles is corresponding to 

Si = AABG, 5*2 = ABCG, 

S^i = ACDG, 5*4 = ADAG, (24) 

as shown in Figure 7, and the total area of this patch 
□ABCD is given by 

S = Si + S2 + S3 + S4. (25) 
Then, patch i receives flux from patch j as much as 



ALij 



aTfSj{nj ■ rji)Si{ni ■ rtj) 



27r\rji\^ 



(26) 



where 

rij = re J - rci, (27) 

rji = rca - rcj, (28) 

as illustrated in Figure 8. We sum up all combinations of 
pair patches on the condition of 

(nj-rji)>0, and (n^ • n^) > 0, 

and obtain the total flux from all patches, i.e.. 



j k 

Here, we set 



2n\rji\^ 



(29) 



(30) 



10, 



(31) 



for not-blocked 
for blocked 

depending whether the line connecting the centers of grav- 
ity is blocked by patch k as shown in Figure 8. 

The surface (patch i) absorbs photon flux and reemits 
by a black-body spectrum with a local temperature of Tj, 
which is determined by 

crT^Si = rjBBLi + aT^^^^^Si, (32) 

where ?7bb is the efficiency of the irradiation, T^^org the 
unheated temperature of the surface. We adopt 50% effi- 
ciency, i.e., 77BB = 0.5 unless otherwise specified. 

Using the window functions, WV(A) and M^b(A), for V- 
and _B-filters (Allen 1973), respectively, we have calcu- 
lated V-, B-, and bolometric magnitudes of the patches. 



I.e., 



AEi{Y) = im ■ e)Si / Wv{X)Bx{Ti)d\, (33) 
Jo 

AEi{B) = {m ■ e)Si / WB{X)Bx{Ti)d\, (34) 
^0 



AEi{ho\) = {m ■ e)SiaTt, 



(35) 



respectively. Here, the unit vector e denotes the direction 
from the patch to the Earth, B\{Ti) the Planck function 
for the temperature of Tj. We have summed up the patches 
having (n^ • e) > 0, i.e., 

E{V) = Y,^Ei(y), for (ni-e)>0 (36) 

i 

E{B) = AEi{B), for {m • e) > (37) 

i 

£(bol) = ^ A£;i(bol), for {m • e) > 0. (38) 
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We have finally calculated the ^-magnitude, my, B- 
magnitude, tub, and the bolometric magnitude, mboh as 

= ^ (- log 33.582) + my,o + 5 log(^), (39) 
"^B = ^(- log £;(B)+ 33.582) + mB,o + 5 log(^), (40) 

mboi = \{- log^(bol) + 33.582)+4.75+5 log(^), (41) 

where d is the distance to the object in units of pc, nivfi 
is a constant determined to satisfy 

mv = mboi, at T = 6500 K, (42) 

and mB,o is also a constant determined to satisfy 

ms = my, at r = 9600 K (AO), (43) 

so that we use 

my,o = 4.52, (44) 

and 

ms.o = 4.91. (45) 

In addition, in order to compare our brightness with the 
recent CCD photometry, we also calculate the Cousins 
Rc- and /c-magnitudes with the window functions given 
by Bessell (1990) along the same way as those described 
above for V- and B-magnitudes. Then, we adopt 

ruR^Q = 4.91, (46) 

and 

m/,0 = 4.50. (47) 

The most time-consuming part of the calculation is pro- 
portional to cubic of the patch number N, that is, oc . 
Therefore, it is hard to increase the patch number sub- 
stantially more than the present maximum case, i.e., 

N = 32 X 64 + 32x64+ (32x64x2 + 8x64) = 8704, (48) 

because it takes about a half day on the Compaq (old 
DEC) Alpha 21264 (500 MHz) machine but about two 
days on the Intel Pentium III (500 MHz) machine to com- 
pute one time sequence for the 1999 U Sco outburst, about 
3000 steps for the my and tub light curve calculations 
(Hachisu et al. 2000a). 

5. T CORONAE BOREALIS 

5.1. Nature of T CrB Outbursts 

T Coronae Borealis (T CrB) is one of the well observed 
recurrent novae and characterized by a secondary maxi- 
mum occurring 150 days after the primary peak. His- 
torically, T CrB bursted twice, in 1866 and 1946, with the 
light curves very similar each other (e.g., Pettit 1946a, b, 
c, d). Large ellipsoidal variations of optical and infrared 
light curves in quiescence suggest that an M3 giant fills its 
Roche lobe (e.g., Bailey 1975; Lines et al. 1988; Yudin & 
Munari 1993; Leibowitz et al. 1997; Shahbaz et al. 1997; 
Belczynski & Mikolajewska 1998). From the ellipsoidal 
variations, the orbital period and ephemeris were deter- 
mined to be HJD 2,431,931.05 + 227.67^; at the epoch of 
spectroscopic conjunction with the M-giant in front (Lines 
et al. 1988). 

There have been debates on the nature of the hot com- 
ponent of the binary system. Radial velocity curves of 
Ki 31 km s~^ {hot component) and K2 ^ 23 km 



s ^ {cool component) indicated the mass ratio of g = 
M2/M1 - 1.4 (Kraft 1958; Paczynski 1965; Kenyon & 
Garcia 1986). Adopting the inclination angle oii = 68°, 
the resulting masses were Mi = Mhot ~ 1.9Mq and 
M2 = Afcooi ^ 2.6Af0, thus placing the hot compo- 
nent above the Chandrasekhar limit. Therefore, Webbink 
(1976) and Webbink et al. (1987) proposed an outburst 
mechanism of T CrB based on their main-sequence accre- 
tor model: the outburst was caused by the transfer of a 
chunk of matter ejected by the giant and. then, by the 
subsequent shock dissipation of the kinetic energy. This 
accretion event model was investigated further by 3D nu- 
merical simulations (Cannizzo & Kenyon 1992; Ruffert, 
Cannizzo, & Kenyon 1993). 

However, Selvelli, Cassatella, & Gilmozzi (1992) have 
been opposed to the main-sequence accretor model from 
their analysis of lUE data in quiescence, which indicates 
the existence of a mass-accreting white dwarf (WD) as fol- 
lows: (1) the bulk of the luminosity is emitted in the UV 
with little contribution to the optical; (2) the presence of 
strong He II and N V emission lines suggests tempera- 
tures of the order of 10^ K; (3) the rotation broadening of 
the high-excitation lines is larger than 1300 km s^^ with- 
out inclination correction. Moreover, the detection of X- 
rays (Cordova, Mason, & Nelson 1981) and the presence 
of flickering in the optical light curves (e.g., lanna 1964; 
Zamanov & Bruch 1998) are also naturally explained in 
terms of accretion onto a WD. 

The mass accretion rate estimated by Selvelli et al. 
(1992) in quiescence is very high (Mace ~ 2.5 x 1O~^M0 
yr""'^) and is exactly required by the thermonuclear run- 
away (TNR) theory to produce a TNR event every 80 yr on 
a massive (> I.3M0) WD. They also pointed out that the 
measurement of the emission-line radial velocities (Kraft 
1958) was intrinsically difficult and uncertain. Thus, the 
1866 and 1946 outbursts can be interpreted in terms of a 
TNR event on a very massive WD (e.g., Starrfield et al. 
1985). 

5.2. Very Massive White Dwarf in T CrB 

Rapid decline rates of the light curves indicate a very 
massive WD close to the Chandrasekhar limit. Assuming 
the solar metallicity {Z = 0.02) of the WD envelope, Kato 
(1995, 1999) calculated nova light curves for the WD mass 
of MwD = 1-2, 1.3, 1.35 and 1.377 Mq, and foimd that the 
light curve of the 1.377M0 model is in better agreement 
with the observational light curve of T CrB than the other 
lower mass models. 

Recently, other observational supports of a massive WD 
in T CrB have been reported. Belczynski & Mikolajewska 
(1998) derived a permitted range of binary parameters, 
MwD = 1-2 ± O.2M0 and q = Mrq/Mwd ^ 0.6, i.e., 
Mrg ^ 0.6 — O.8M0, from amplitude of the ellipsoidal 
variability and constraints from the orbital solution of 
M-giants. In Shahbaz et al. (1997), a massive WD 
of MwD = 1.3 — 2.5M0 is suggested from the infrared 
light curve fitting. Moreover, the measurement of the 
radial velocities (Kraft 1958) has been revised includ- 
ing the effect of gas-streams around the hot component 
(Hric et al. 1998). Hric et al. (1998) estimated the 
masses, Mwd = 1-2 ± 0.2 Mq for the hot component and 
-^RG = 1.38 + 0.2 Mq for the M-giant, assuming the in- 
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clination of i = 68°. Combining these permitted ranges of 
the WD mass in T CrB, we may conehide that a mass of 
the WD is between Mwd = 1.3— 1.4M0, which is very con- 
sistent with the hght curve analysis Mwd ~ 1.37— 1.38M0 
by Kato (1999). 

5.3. Secondary Maximum of Outbursts 

The secondary maximum in outbursts is rarely observed 
in fast novae and in other recurrent novae. Selvelli et al. 
(1992) suggested a possibility of an irradiated stationary 
shell around the binary system, although the presence of 
such a shell is just speculation. Recently, Hachisu & Kato 
(1999) proposed another mechanism of the second peak, 
i.e., radiation-induced tilting disk around a massive WD. 
The main results of their analysis are: (1) the first peak 
is naturally reproduced by a fast developing photosphere 
of the WD envelope based on the TNR model of a very 
massive WD, Mwd ~ 1.35M0; (2) the second peak is not 
fully reproduced by an irradiated M-giant model as sim- 
ply estimated by Webbink et al. (1987); (3) instead, the 
second peak can be well reproduced if an irradiated tilting 
accretion disk around the WD is introduced together with 
the partly irradiated M-giant companion. 

Such a radiation-induced, tilting instability of an accre- 
tion disk has been suggested by Pringle (1996) when a 
central star is as luminous as the Eddington limit. This 
radiation-induced instability sets in if the condition 



M«., 



10-7 M© yr-i 



< 



-Rdisk 
10 Rq 

i?WD 

0.003 Rq 



1/2 



Lhoi 



2 X 10^8 ergs 



1/2 



M 



W^D 



1.37 MfT 



-1/2 



(49) 



is satisfied (Southwell, Livio, & Pringle 1997), where i?disk 
is the radius at the edge of the optically thick accretion 
disk. Selvelli et al. (1992) estimated the accretion rate of 
T CrB as low as Mace ~ 0.25 x IQ-'^Mq yr'^ which meets 
condition (49) above even if the optically thick radius is 
as small as iZdisk ~ 1 Rq- Therefore, it is likely that the 
radiation induced instability grows during the outburst in 
T CrB. The growth timescale of warping is estimated by 
Pringle (1996) and Livio & Pringle (1996) as 



~ 30 



M, 



disk 



10-8 Mq 

-^bol 



M 



WD 



1.37 MfT 



-Rdisk 



R. 



-1/2 



2 X 10^^ ergs s" 



day 



(50) 



where we assume that the a-parameter in the Shakura- 
Sunyaev (1973) standard accretion disk is a ~ 0.1. Thus, 
the growth timescale is short enough to excite warping of 
the accretion disk. 

5.4. Revised Model of Warping Accretion Disk 

We have revised Hachisu & Kato (1999) model to be 
consistent with the observations described below. The 
main difference is in the accretion disk, which plays an 
essential role in the mid/latc phase of the V light curve. 
Hachisu & Kato (1999) assumed that the size of the ac- 
cretion disk is as large as 0.7 times the Roche lobe size, 
i.e., iidisk ~ 50 Rq. This is because the mass transfer is 
not a wind- fed type but a Roche lobe overflow type (e.g., 



Schandl et al. 1997). However, SelvcUi et al. (1992) es- 
timated the optically thick region of the accretion disk in 
quiescence is as small as Rdisk ^ 1 Rq because the disk lu- 
minosity contributes mostly to the satellite UV. Belczyhski 
& Mikolajcwska (1998) pointed out that the optical con- 
tribution of the accretion disk is clearly visible in the U, B, 
and V bands in 1981-85 (optically high state) while it is 
practically absent in 1990-94 (optically low state) as first 
shown by Payne-Gaposchkin & Wright (1946). Thus, Bel- 
czyhski & Mikolajcwska derived the disk radius (hot spot) 
of i?disk ~ 0.1a (20 Rq) in the optically high state. 

Therefore, we here adopt a small size of the optically 
thick disk, i.e., i?disk ^ 1 Rq just before the outburst. A 
large accretion disk of ^ 50 Rq requires that (1) the ir- 
radiation efficiency is too small and that (2) the disk is 
inclined from too early phase of the outburst. If we intro- 
duce a small size of the accretion disk, we do not need to 
assume such unnatural assumptions. This accretion disk 
satisfies the warping condition (49) so that we adopt an 
analytic form proposed by Maloney et al. (1996). 

5.5. Light Curve of Early Phase 

The very fast decay during the first 10 days can be well 
reproduced only by a WD photosphere of a very massive 
WD, Mwd ~ 1.35—1.37 Mq, as already shown by Kato 
(1999) and Hachisu & Kato (1999). Wc have found that 
the decay rate of the 1.37 M© light curve is in better agree- 
ment with the observational points than the others. The 
difference from Hachisu & Kato (1999) comes from the 
difference of the disk size which affects the light curve 10 — 
30 days after maximum. Hachisu & Kato (1999) adopted 
a large accretion disk which in general contributes more 
to the luminosity. Therefore, they reduced the WD lumi- 
nosity to be consistent with the observational points by 
occulting the WD photosphere with the accretion disk. In 
the present model, on the other hand, the small disk can- 
not occult a large part of the WD photosphere during the 
day 10 — 30. Thus, we need a more massive WD to gen- 
erate a faster decline. Therefore, we have determined the 
WD mass of Mwd = 1-37 ± 0.01 M© in T CrB, giving a 
much better fit with the observational points even in the 
day 0—10. 

Deutsch (1948) suggested from the observations of Ash- 
brook (1946), Wright (1946), and Pcttit (1946a) that the 
color index of T CrB was about zero February 12 (4 days 
after maximum). Our model indicates the photospheric 
color of {B-V)o = -0.2 (Tph ~ 18,500 K, iiph ~ 18 Rq), 
which is consistent with the observational color because 
the reddening is about E{B — V) = 0.15. 

The mass of the companion is estimated as follows: 
Adopting Mwd = 1-37 M© and the radial velocity of the 
cool component K2 = 23.32 km (Kenyon & Garcia 
1986), i.e.. 



f{M) 



= 0.30±0.01M, 



0: 



(51) 



(Mwd + Mrg) 
we have Mrg = 1-30, 1.25, 1.16, and 1.0 M©, for the 
inclination angle of i = 70°, 68°, 65°, and 60°, respec- 
tively. Following Belczynski & Mikolajcwska (1998), 
q = Mrg /Mwd ^ 0.6 and an inclination angle of z ^ 60°, 
we have assumed the red giant mass of Mug = 1-0 M© {i = 
60°). So we show here only the results for Mrg = l.OM©. 
In this case, the separation is a = 209. 2i?©, the effective 
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radius of the inner critical Roche lobe for the WD compo- 
nent is Rl = 85.Oi?0, the effective radius for the red giant 
companion is R2 = R2 = 73.6-R0. We have examined the 
other two cases of the red giant mass, i.e., Mrg = 1-25 Mq 
(i = 68°) and A/rg = 1.16 Mq [i = 65°), and obtained 
similar light curves to that of Mrq = 1.0 Mq {i = 60°). 

Figure 9 shows the best fitted model to the 1946 out- 
burst together with the observational points. Here, we 
assume the solar metallicity {Z = 0.02) and hydrogen 
content oi X ~ 0.70 for the WD envelope. We also as- 
sume that the V magnitude reached its maximum on HJD 
2,431,861. The dotted curve denotes the total V light 
from the WD photosphere and the non-irradiated red gi- 
ant at TRG.org — 3100 K. Ellipsoidal light variations can 
be seen in the late stage (dotted). The effective tempera- 
ture of the M-giant is estimated to be TRG,org = 3100 K 
from our light-curve fitting, which is a little bit lower than 
Tcff = 3560 K by Belczyhski & Mikolajewska (1998), but 
roughly consistent with Tcff' = 3200 K by Shahbaz et al. 
(1999). 

If we include irradiation of the companion, the light 
curve changes to the dash-dotted curve. This irradiation 
effect can be well seen both 10 — 40 days after the optical 
maximum and 120 260 days after the optical maximum. 
The irradiation effect was once examined and discarded by 
Webbink et al. (1987) as a main mechanism of the sec- 
ondary maximum, mainly because its maximum largely 
deviates from the observational secondary maximum as 
seen in the figure. 

The solid curve includes further the irradiated accre- 
tion disk. Here, we have started the calculation by as- 
suming the disk size of a = 0.007 in equation (11), i.e., 
.Rdisk = 0.6 Rq. 

The photospheric radius expands to iiph ~ 100 Rq in 
order to fit the peak luminosity of the V light. We have 
derived the distance modulus of (m — M)v = 10.2. The 
light curve during the first 10 days of the 1946 outburst 
can be well reproduced only with the fast-developing WD 
photosphere of Mwd = 1-37 Mq as already shown by 
Kato (1999). The envelope mass of the WD at the op- 
tical maximum is AM = 3.3 x IO^^A/q from our wind 
solutions. Then, the mass accretion rate was estimated to 
be Mace = 4 X 1O~^M0 yr~^ during the quiescent phase 
between 1866 and 1946. This value is roughly consistent 
with the estimation of 2.5 x 10~^Mq yr^^ by Selvelli et 
al. (1992). 

After the WD photosphere shrinks to smaller than 
.Rdisk ~ 0.6 Rq, the accretion disk reappears but its contri- 
bution to the V light is small as shown in Figure 9, because 
radiation emits mainly in UV. We assume a flaring-up rim 
of the accretion disk, i.e., (3 = 0.30 in equation (12), dur- 
ing the strong wind phase for the 1946 outburst of T CrB, 
since Hachisu et al. (2000a) suggested that the rim of the 
accretion disk is flaring-up during the 1999 outburst of U 
Sco, Then, the V light is reduced 20 — 30 days after maxi- 
mum because the disk rim makes shadow on the red giant 
component. The warping of the accretion disk does not 
occur yet, because the strong wind suppresses the warp- 
ing. The spherically symmetric wind momentum is much 
larger than the radiation momentum as shown by Kato & 
Hachisu (1994). 



5.6. Light Curve Model of the Second Peak 

Warping (or wobbling) instability of the accretion disk 
begins to grow after the strong wind stops ^ 60 days af- 
ter maximum (on HJD 2,431,920) and, as a result, the V 
light attains another peak. We expect that the radiation 
induced warping of the accretion disk grows in a time scale 
of a few tens of days as derived in equation (50) and further 
that the mass transfer from the red giant to the WD in- 
creases because the red giant surface is strongly irradiated 
and the Roche lobe overflow is enhanced. Such enhance- 
ment of the mass transfer certainly increases the radius of 
the optically thick region of the accretion disk up to ~ 10 — 
20 Rq as suggested by Belczynski & Mikolajewska (1998). 
Thus, we conclude that the radius of the optically thick ac- 
cretion disk increased from a = 0.007 (i?disk ~ 0.6 Rq) to 
a = 0.07 (i?disk ~ 6 Rq) during the day -60 - 150, by 
fitting with the observational points as shown in Figure 9. 

The second peaks of the 1866 and 1946 outbursts are 
very similar. This is explained as follows: (1) The first 
peak ends up when the wind stops about 60 days after 
maximum. This period (— 60 days), between the opti- 
cal peak and the termination of the strong wind, is al- 
most uniquely determined by the WD mass. (2) The sec- 
ond peak begins to rise when the warping of the accretion 
disk grows, which is about 120 days after maximum. The 
growth time is estimated to be about a few times 30 days 
as in equation (50). This growth time depends on the mass 
(Mdisk) and the optically thick radius (i?disk) of disk, which 
are determined mainly by the mass transfer rate from the 
red giant companion. We can naturally expect the similar 
time scales between the 1866 and 1946 outbursts, because 
the irradiation of the red giant by the WD photosphere is 
almost the same. (3) Thus, we are able to obtain quite 
similar light curves if the orbital phases of the 1866 and 
1946 outbursts happen to be the same (see the light curves 
of the 1866 and 1946 outbursts, Pettit 1946d). Acciden- 
tally, the two orbital phases were almost similar to each 
other. 

The precessing period of the warping disk is similar to 
the growth time of warping (e.g., Livio & Pringle 1996). 
About 60% faster precessing angular velocity than the or- 
bital motion is required from the phase relation between 
the rising shoulder of the second peak near HJD 2,431,960, 
a small dip near HJD 2,432,090, and then a small bump 
near HJD 2,432,120 as shown in Figure 9. These dips are 
caused by a large shadow on the companion cast by the 
accretion disk. In other words, the accretion disk blocks 
the light from the WD photosphere and makes its shadow 
on the companion surface. Finally, we have determined 
-Pprec = 145 days by fitting. This period is similar to the 
164 days precessing period of SS 433. 

Webbink et al. (1987) summarized the observational 
development of the 1946 outburst, showing the color in- 
dex at the second peak of B — V ^ 0.4 after correcting the 
reddening of E{B — V) = 0.15. They further pointed out 
that the refiection by the M-giant is about B — V ~ 1.35, 
therefore, rejected the possibility of reflection by the M- 
giant. Instead, they suggested that the much bluer color 
index of ~ 0.4 is typical of the accretion disks seen in out- 
bursting accretion-powered symbiotic stars (e.g., Kenyon 
& Webbink 1984). In our warping-disk model, the color 
index is ranging from B-V = 0.15 (on HJD 2,432,000) 
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to B~V = 0.35 (on HJD 2,432,050), which is reasonably 
consistent with the observational color at the second peak 
because the reflection effect by the warping disk dominates 
during the second peak rather than the irradiated M-giant. 

Ever since the emission lines of O III, N III, He II, [O III], 
[Ne III], and [Ne V] first appeared, they have been always 
present until two years after the outburst, except in early 
June, at the most, from May 15 to July 17, 1946 (~ 105— 
135 days after maximum, Sanford 1949). In our model- 
ing, the warping-up disk occulted both the WD and the 
innermost region of the disk during these period, being 
consistent with the disappearance of the emission lines. 

5.7. The Distance, Envelope Mass of the White Dwarf, 
and Possibility of SN la Explosion 

The distance to T CrB is estimated to be d = 0.94 
kpc with the absorption of Ay = 0.35 (Harrison, John- 
son, & Spyromilio 1993) because the distance modulus 
is (to - M)o = (to - M)v - Ay = 10.2 - 0.35 = 9.85. 
This is consistent with both Harrison et al.'s estimation of 
d = 1100 pc, Belczyiiski & Mikolajewska's (1998) estima- 
tion of d = 960 ± 150 pc, and Bailey's (1981) estimation of 
d = 1180 pc. Harrison et al. determined the infrared spec- 
tral type of the cool companion as M3 III. Assuming the 
absolute magnitude of Mk ~ —5.3 (M3 III), they obtain 
the distance modulus of (m — M)o — rriK — Mk — Ak = 
4.79 - (-5.3) - 0.04 = 10.05, where they used a relation 
between the absorptions of Ak and Ay given by Rieke & 
Lebofsky (1985), i.e., AK/Ay = 0.112. 

The envelope mass of the WD at the first peak is 
AM = 3.3 X 1O~^M0 from our wind solutions. There- 
fore, the mass accretion rate can be estimated to be 
Ma_cc = 4 X 10~^A'/q yr~^ if the envelope mass had been 
accumulated during the quiescent phase between 1866 and 
1946. This value is roughly consistent with the estimation 
of 2.5 X IO-^Mq yr-i by SelveUi et al. (1992). The wind 
mass loss rate; at the first peak is M„ind = 3.7 x 1O~^M0 
yr"^ while the nuclear burning rate is as small as Mnuc = 
9.3 X 1O~^M0 yr~^. The wind mass loss rate is quickly de- 
creasing to Afwind ^ 1 X yr~^ about 30 days after 
the first peak. Then, 90% of the envelope mass has been 
blown in the wind. The residual 10% is left on the WD as 
a part of the helium layer. Therefore, the net growth rate 
of the WD is Mhc 4 x IO^'U'/q yr^^. 

At the end of this section, we discuss whether or not 
T CrB will explode as a Type la supernova in the near 
future. The conditions for SN la explosions are summa- 
rized as follows (see, e.g., Nomoto & Kondo 1991): (1) 
the accretion rate of the C-l-0 WD should be higher than 
the critical accretion rate under which helium detonation 
occurs, i.e.. 

Mace > Mne-det ~ 1 X IQ-^ Mq yr'^, (52) 

for Population I stars; (2) the initial mass of the C-l-0 WD 
at the beginning of the mass transfer is less massive than 

Mc+ofl < 1.2 Mq. (53) 
T CrB satisfies the first condition as has already been sug- 
gested to be Mace 4 X IO^^Mq yr'^ (even if Mace ~ 
2.5 X 10-^Mq yr-i by SelveUi et al. 1992). For the sec- 
ond condition, we cannot estimate the initial WD mass of 
T CrB. However, evolutionary scenario strongly suggests 



that the mass of the WD in T CrB has been increased up 
to the present mass from a much lower initial mass (see 
e.g., the evolutionary path to SNe la proposed by Hachisu 
et al. 1999a). Moreover, it is theoretically predicted that 
an upper limit of degenerate carbon-oxygen core is about 
1.07 Mq, i.e., Mc+o,o < 1-07 M© (Umeda et al. 1999b). 
Thus, we may conclude that T CrB will explode as a Type 
la supernova in quite a near future if the hot component 
is a carbon-oxygen white dwarf. 

5.8. Summary of T CrB Outburst 

The main features of the T CrB outbursts can be well 
understood with a thermonuclear runaway model on a very 
massive white dwarf having a warping accretion disk. We 
have summarized the following important points from the 
fitting of the modeled light curve with observation: 

(1) The early phase (^ 10 days) visual light curve of T 
CrB outbursts can be well reproduced only by a bloated 
white dwarf photosphere of Mv^^d = 1.37 ± 0.01 Mq for 
solar metallicity oi Z = 0.02. This is because decline rates 
of optical light curves depends sensitively on white dwarf 
mass especially for white dwarf masses near the Chan- 
drasekhar mass limit. Thus, we are able to specify the 
mass of white dwarf with a great accuracy (three digits) . 

(2) The ensuing early phase (~ 10 — 40 days) visual light 
curve cannot be reproduced by any combination of a white 
dwarf photosphere and a non-irradiated red giant. It re- 
quires an irradiation effect of the red giant component. 
This heated-up red giant surface is consistent with the 
very early appearance of M-type spectrum (TiO bands) 
shortly after the outburst, indicating the luminosity of M- 
type giant much brighter than 10th magnitude (Deutsch 
1948). 

(3) Through the mid and late phase (^ 40 300 days), the 
secondary maximum can be well modeled if we introduce 
the radiation-induced warping instability of the accretion 
disk, which sets in after the wind stops ~ 60 days after 
maximum. Introducing an analytic form of warping (Mal- 
oney et al. 1996), we obtain the size of the accretion disk 
as large as i?disk ~ 6 Rq and the precessing period of the 
warping disk as fast as ~ 140 days to reproduce the late 
phase light curve. 

(4) The envelope mass of the bloated white dwarf is 
3.3 X 1O~^M0, suggesting the mass accretion rate of 
0.4 X 10-^Mq yr"^ between the 1866 and 1946 outbursts. 

(5) The wind carries away about 90% of the envelope mass, 
so that the residual 10% can be accumulated to the helium 
layer on the white dwarf. As a result, the white dwarf 
grows at a rate of 0.4 x 10~*M© yr~^. 

(6) T CrB will certainly explode as a Type la supernova 
in quite a near future if the white dwarf consists of carbon 
and oxygen. 

6. RS OPHIUCHI 

RS Oph is also one of the well-observed recurrent novae. 
It is characterized by a long orbital period of 460 days (Do- 
brzycka & Kenyon 1994, but recently 457 days by Fekel et 
al. 2000), a relatively short recurrence period of ~ 10 — 20 
yrs compared with 80 yrs of T CrB. It has been suggested 
that a companion (M-giant) star is underfilling the Roche 
lobe and losing its mass by massive stellar winds (e.g., 
Dobrzycka et al. 1996). RS Oph underwent five recorded 
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outbursts (in 1898, 1933, 1958, 1967, and 1985), with the 
hght curves being very similar to each other (e.g., Rosino 
1987). The latest (1985) outburst has been observed at 
all wave lengths from radio to X-rays (sec papers in Bode 
1987). 

Although there had been intense debates on the mecha- 
nism of RS Oph outbursts (e.g., Livio, Truran, & Webbink 
1986; Webbink et al. 1987), various observational as- 
pects favor TNR models on a very massive WD (see, e.g., 
Anupama & Mikolajcwska 1999, for a recent summary). 
Rapid decline rates of the V light curves also indicate a 
very massive WD close to the Chandrasekhar limit. Kato 
(1991) has first calculated RS Oph outburst light curves 
for the WD masses of 1.33, 1.35, 1.36 and 1.37 Mq and 
found that the light curve of the 1.36AfQ model is in better 
agreement with the observational light curve of RS Oph 
than the other mass models. 

Since Kato (1991) adopted an old opacity in her calcula- 
tion, Hachisu & Kato (2000a) recalculated the light curves 
of RS Oph outbursts with the new opacity (e.g., Iglesias 
& Rogers 1996), including irradiation effects of both an 
M-giant companion and an accretion disk around the WD. 
Their main results are summarized as: (1) the first 4 days 
light curve of the outbursts can be well reproduced with 
a TNR model based on a Mwd = 1-35 ± O.OIM© and an 
apparent distance modulus of (m — M)y = 11.09. (2) 
To reproduce the V and UV light curves in the late stage 
{t ~ 4 100 (lays after the optical maximum), the contri- 
butions of the irradiated M-giant and the irradiated ac- 
cretion disk must be introduced. (3) The distance to RS 
Oph is estimated to be 0.6 kpc both from the V and UV 
light curve fittings with the absorption of Ay — 2.3. (4) 
About 90% envelope mass of the WD has been blown off in 
the wind but the residual 10% has been left on the WD. 
Thus, Hachisu fc Kato (2000a) has concluded that RS 
Oph will explode as an SN la in the near future if the hot 
component is a carbon-oxygen white dwarf. 

In Hachisu & Kato's (2000a) calculations, they assumed 
the solar metallicity {Z = 0.02) of the M-giant enve- 
lope. However, it has recently been suggested that the 
cool component of RS Oph is a metal-poor giant. Scott ct 
al. (1994) obtained the carbon abundance of [C/H]= —3 
for RS Oph 143 days after the 1985 outburst and of 
[C/H]= —4 in 1992. Similar depletions of carbon, oxy- 
gen, neon, sulphur, and iron with respect to solar values 
arc suggested by Contini, Orio, & Prialnik (1995). They 
obtained one or two tenths of the solar values for C, O, 
Ne, S, and Fe, in order to reproduce the line ratios based 
on their shock models of nova shells. Smith et al. (1996) 
estimated metal abundance of AG Dra as [Fe/H]= —1.3 
and suggested that all yellow symbiotic stars including RS 
Oph are metal-poor. They attributed K6 I-II to the cool 
component of RS Oph. 

If the envelope of the white dwarf is also metal-poor, the 
mass determination of the white dwarf from the light curve 
fitting may not be correct because the optically thick wind 
is driven by iron lines and the decay timescale is much 
longer in low metallicity envelopes. Therefore, we must 
reexamine the light curves of the 1985 outburst by intro- 
ducing lower metallicities and study the possible range of 
the metallicity and the WD mass. 



6.1. The Model of the 1985 Outburst 

Assuming an extremely massive white dwarf of Mwd = 
1.377 Mq just prior to the SN la explosion, we have cal- 
culated light curves for various metallicities oi Z = 0.001, 
0.004, 0.01, and 0.02 as shown in Figure 10. The V light of 
the first five days (HJD 2,446,092 - 2,446,096) can be re- 
produced only by the bloated WD photosphere as already 
shown in Hachisu & Kato (2000a). The case oi Z = 0.001 
is too slow to be compatible with the observations. Thus, 
we may conclude that the metallicity is Z > 0.004. Here, 
we assume Z = 0.004 for RS Oph as a lower limit of metal- 
Hcity. For the set of Mwd = 1-377 Mq and Z = 0.004, we 
have obtained the envelope mass at the optical maximum, 
AM = 2.2 X 1O~^M0, which indicates a mass accretion 
rate of 1.2 x 1O~'^M0 yr~^ during the quiescent phase be- 
tween the 1968 and the 1985 outbursts. 

The irradiation effects of the accretion disk and the red 
giant (RG) are necessary to reproduce the optical and UV 
light curves of the outburst as already shown by Hachisu 
& Kato (2000a). Dobrzycka & Kenyon (1994) obtained 
the radial velocities of the M-typc feature and the A-type 
feature but concluded that the A-type absorption feature 
does not associate with any of the components. Thus, RS 
Oph is a single line binary with 
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(54) 

(Mwd + Mrg)" 
from the semiamplitude of the radial velocity Krq = 
12.8 ± 2.0 km s"^ and the orbital period of P = 460 
days. They suggested a relatively low inclination angle 
of'?; - 30—40°. If Mwd = 1-377 Mq, then we obtain 
Mrg = 0.5 Mq for i = 30°, Mrg = 0.7 Mq for i = 33°, 
Mrg = 0.9 Mq for i = 36°, Mrg = 1.0 Mq for i = 37°, or 
Mrg = 1.2 Mq for i = 39°. We assume Mrg = 1.0 Mq 
and i = 37° in the present paper, although we have al- 
ready confirmed that the other set of Mrg = 0.7 Mq and 
i = 33° reproduces almost the same light curves. 

For the cool component of the RS Oph system, as sug- 
gested by Dobrzycka ct al. (1996), the RG lies well within 
the inner critical Roche lobe, i.e., 

Rrg=iR; (7<1), (55) 
where i?2 is the effective radius of the inner critical Roche 
lobe for the RG component and 7 is a numerical factor. 
Dobrzycka et al. (1996) adopted 7 = 0.4 for the distance 
oi d = 1.5 kpc to RS Oph. On the other hand, Hachisu & 
Kato (2000a) estimated the distance of d = 0.6 kpc to RS 
Oph much shorter than Dobrzycka et al.'s value, so that 
we adopt 7 = 0.25 for d = 0.6 kpc in this paper. Even 
if we use 7 = 0.4, the light curve itself is not so differ- 
ent from that for 7 = 0.25 as shown in Hachisu & Kato 
(2000a). Only the difference is in the irradiation effect of 
the M-giant, but the irradiation itself is relatively small 
compared with the case of T CrB (7 = 1). 

Recent studies on the radii of M-giants in symbiotic stars 
indicate 7 ^ 0.5 for the normal S-type symbiotic stars 
(e.g., Miirset & Schmid 1999). The cool component of 
the RS Oph system is not a normal giant. It is suggested 
to be metal-poor and sometimes classified to K-giant (e.g.. 
Smith et al. 1996) rather than M-giant. Therefore, we 
may adopt the relatively small radius of 7 = 0.25. 

As for the orbit of the binary system, the ephemeris 
for the inferior conjunction of the red giant in front is 
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2,444,999.9+460 x E (Dobrzycka & Kenyon 1994). The 
separation is a = 325. Oi?©, the effective radii of the inner 
critical Roche lobes for the WD component and the RG 
component are Rl = 138.6i?0 and = lO8.2i?0. For 
7 = 0.25, we have Rrg = 0.25i?^ ~ 27Rq. 

We assume 50% efficiency of the irradiation of the red 
giant (RG) component (t/rg = 0.5). The nonirradiated 
photospheric temperature of the RG is a parameter which 
is determined to be fitted with observation. Here, we adopt 
^RG = 3400 K against 7 = 0.25 as shown in Figure 11. 

We have also included the contribution of irradiated ac- 
cretion disk (ACDK) luminosity. Here, we also assume 
50% efficiency of the ACDK irradiation (77DK = 0.5). The 
viscous heating is neglected because it is much smaller 
than that of the irradiation effects. The temperature of 
the unheated surface of the ACDK including the rim is 
assumed to be Tdisk = 2000 K. We have checked two other 
cases of Tdisk ~ 1000 and K and have found no significant 
differences in the light curves. The power of v is assumed 
to be v = 9/8 from the standard disk model. We have 
checked the dependency of the light curves on the param- 
eter v by changing from = 9/8 to 1/ = 2 but cannot find 
any significant differences as far as the disk rim is small 
(/?'= 0.01—0.05). 

The late-phase V light curve {t ^ 5 — 100 days after 
maximum) indicates strong irradiation of the ACDK. We 
have examined the case of a lobe-filling RG companion 
(7 = 1), but the irradiation of the RG is too luminous 
to be compatible with the observational light curves. On 
the other hand, a smaller size of the RG such as 7 = 0.25 
(^ 27 Rq) gives a reasonable fit with the observations as 
shown in Figures 11 and 12. 

To fit the late-phase light curve, we have calculated the 
two cases of the disk parameters, i.e., {a, (3) = (0.1, 0.05) 
and (0.008, 0.05), as shown in Figure 11. Here, two solid 
lines denote the cases of a = 0.1 (i?disk = 14 Rq) and 
(3 = 0.05 and a = 0.008 (i?disk = 1 Rq) and /3 = 0.05, 
respectively. Something between these two can roughly 
reproduce the light curve of the 1985 outburst but the 
gradual transition from a = 0.1 to a = 0.008 is much bet- 
ter for reproducing the observations as shown in Figure 
12. It is very likely that the surface of the accretion disk 
has been dragged by the strong wind and been gradually 
blown in the strong wind. As a result, the accretion disk 
has become smaller during the strong wind phase (~ 70 
days after maximum). Thus, we adopt 




(t-to)/72 

, ioTto<t<to + 72 day (56) 



where ao = 0.1, ai = 0.008, t is the days after maximum, 
to = 7 days is the reappearance day of the accretion disk, 
and 72 days is roughly the wind period. 

We have also calculated UV light curves with a response 
function of 911 — 3250A to fit the reddening corrected UV 
data (Snijders 1987a). Figure 11 shows two cases of UV 
light curves with the parameters of (a = 0.1, (3 = 0.05) 
and {a = 0.008, /3 = 0.05). The calculated UV light curves 
can reproduce well the UV observations both for a = 0.1 
and a ~ 0.008 with the distance of d = 0.57 kpc because 
the innermost region (< 0.5 Rq) of the accretion disk ra- 
diates strong UV but the outer region of > 1 Rq mainly 
contributes to optical light. 



6.2. Development of the 1985 Outburst 

In this subsection, we summarize the observational de- 
velopments of the RS Oph 1985 outburst and discuss them 
in relation to our model calculation. 

Rosino (1987) and Rosino & lijima (1987) divided the 
outburst into several phases along with the development 
of the optical spectrum: [[Phase — 2 days]] after the opti- 
cal maximum, broad emission lines of H I, He I, and Fe II 
are observed with narrow absorption component (40 — 60 
km s"^) and broad P Cygni troughs (-2700 to -3700 km 
s-i). [[Phase 2 -14 days]], He II 4686 A , [N II], [O I], 
[O III], [Ne III], [S II], and [Fe II] 6374 A appear. [[Phase 
14—29 days]]. He II 4686 A and N III 4640 A become 
strong. [[Phase 29 51 days]], coronal lines appear and be- 
come strong. [[Phase 51 — 72 days]], the coronal line phase 
lasts. [[Phase 72 — 119 days]], the coronal lines disappear. 
Most of high excitation lines, including He II 4686 A, are 
fading. 

Snijders (1987a) presented an overall development of 
the UV light curve during the 1985 outburst of RS Oph as 
already shown in Figures 11 and 12. 

Since RS Oph exploded in a circumstellar envelope 
that is formed by the slow wind of the red giant com- 
ponent (Wallerstein 1958), the rapidly expanding ejecta 
(~ 3,000—4,000 km s-\ e.g., Rosino 1987; Snijders 
1987a) of the outburst is decelerated by the slowly expand- 
ing (< 20 km s~^) circumstellar matter to form a strong 
shock as noticed by Pottasch (1967). The secular varia- 
tion of coronal lines was interpreted by Gorbatskii (1972) 
in terms of a decelerating shock wave (see also Wallerstein 
& Garnavich 1986, for the 1985 outburst). This shock- 
heated layer is also responsible for soft X-ray (Mason ct 
al. 1987) and radio (Hjellming et al. 1986) emission as 
will be mentioned below. 

In our numerical modeling, the photospheric temper- 
ature of the WD envelope has once decreased down to 
Tph = 12, 200 K at the optical maximum (HJD 2,446,093). 
Then, it increases gradually to [[Phase 0—2]] Tph = 12, 200 

— 17,700 K for the WD photospheric radius of Rp\, = 45 

— 21 Rq, [[Phase 2—14]] Tph = 17,700 — 67,100 K 
for i?ph = 21 1.6 Rq, [[Phase 14—29]] Tph = 67,100 

— 114,000 K for i?ph = 1.6 — 0.56 Rq, [[Phase 29— 
51]] Tph = 114,000 — 181,000 K for i?ph = 0.56 — 
0.23 Rq, [[Phase 51—72]] Tph = 181,000 — 302,000K for 
i?ph = 0.23 0.083 Rq, [[Phase 72 119]] Tph = 302. 000 

— 1,020,000 K for i?ph = 0.083 — 0.0037 Rq. The accre- 
tion disk reappears 4 days after the optical maximum, i.e., 
the photosphere shrinks to iiph ^ 8 Rq, which is smaller 
than the size of the accretion disk iZdisk ~ 10 Rq that we 
assume. The optically thick strong wind, which is blowing 
from the WD, stops 77 days after the optical maximum 
(HJD 2,446,170). The photospheric radius shrinks from 
^ph = 0.06 Rq to i?ph = 0.01 Rq within a few days. Cor- 
respondingly, the photospheric temperature increases from 
Tph = 343, 000 K to Tph = 1 X lO*' K. This sharp decrease 
in the photospheric radius causes the sharp drop in the UV 
light curve in Figures 11 and 12, because the effective an- 
gle of cos9 in equation (14) against the disk surface drops 
by a factor of three. The steady hydrogen shell-burning 
also ends 109 days after maximum (HJD 2,446,202). 

The appearance of He II 4686 A , [N II], [O I], [O III], 
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[Ne III], [S II], and [Fe II] 6374 A in [[Phase 2—14]] and 
[[Phase 14 — 29]] is consistent with our model because the 
WD photospheric temperature increases to high enough to 
excite these elements. The appearance of the coronal lines 
are related with the appearance of X-rays heated by the 
blast shock as shown below (Mason et al. 1987) because 
they need high temperature photons of T > 1 x 10^ K as in 
the solar corona. The disappearance of the coronal lines 
is usually understood by the shock breakout at the day 
~ 70. However, the quick disappearance of the coronal 
lines in [[Phase 72 — 119]] might be caused by the extinc- 
tion of the hydrogen shell burning, as suggested by Shore 
et al. (1996), because it is consistent with the extinction 
time of the steady hydrogen shell burning in our model. 

Evans et al. (1988) presented infrared observation (J, 
H, K, and L bands) of the RS Oph 1985 outburst. The 
decline rates of J, H, K, and L are almost the same as 
t2 ^ 30 days (1,2 is the time taken to drop 2 mag from max- 
imum). At each wave length it is noticeable that the IR 
flux of RS Oph had returned to its pre-outburst value by 
the day ~ 85 (after the optical maximum, but Evans et al. 
took the optical maximum on HJD 2,466,095.5 as a origin 
of time) and that the decline continued beyond this level. 
However, the IR fluxes had reverted to their pre-outburst 
levels until the day ^400. The post-outburst minimum 
lasted during about 300 days, i.e., from the day ~100 to 
the day ^400. This post-outburst minimum is discussed 
below in relation to the super-soft X-ray detection on the 
day 251 and the distance to RS Oph. 

Evans et al. (1988) also estimated the IR flux from the 
circumstellar material in the RS Oph system by removing 
the contribution from the M-giant companion. Assuming 
that the He II 1640 A line and the IR continuum do in 
fact arise in the same region, they fitted the IR continuum 
on HJD 2,466,119 by a nebular continuum having electron 
temperature To ^ 1.6 x 10'"' K. In our model calculation, 
we have the photospheric temperature of the bloated WD 
envelope Tph ~ 1.1 x 10^ K on the same day, being roughly 
consistent with Evans et al.'s estimation, though Evans et 
al. attributed it to the cooled shocked layer. Evans et 
al. also estimated that the spectral type of RS Oph cool 
component is K8 HI (±2 sub-divisions) from [2.35] — [2.2] 
color. 

Shore et al. (1996) also divided the outburst into two 
phases along with the development of the UV spectrum: 
The first phase, [[Phase 60 days]] after the optical max- 
imum, is the broad line phase. The high velocity material 
dominated the emission. This is roughly consistent with 
the strong wind phase of our numerical modeling, i.e., until 
HJD 2,446,170. In the second phase, i.e., [[Phase 60—130 
days]], the newly ionized red giant wind dominated the 
spectrum. This period is also consistent with the steady 
hydrogen burning phase of our numerical modeling, i.e., 
until HJD 2,446,202. 

Shore et al. also discussed that the line ratio of Si III] 
1892A to C III] 1910 A, which measures the electron 
density (Nussbaumer & Stencel 1989), steadily decreased 
from an initial value of 1.0 on HJD 2,446,106 to a mini- 
mum of 0.3 on HJD 2,446,164. Then, the ratio began to 
rise thereafter, reaching 0.6 by HJD 2,446,202. This trend 
can be explained by our wind model: (1) The wind mass 
loss rate is gradually decreasing from 10~* to 1O~^M0 



yr ^ and this causes the decrease in the electron density 
because of p ^ Mjvind/47rr^Uwind- (2) After the hot wind 
from the WD stops, the cool red giant wind resumes to 
expand slowly and the newly ionized red giant wind dom- 
inates the spectrum. 

RS Oph became a strong X-ray source with characteris- 
tic temperature of a few times 10^ K (Mason et al. 1987). 
Mason et al. observed X-rays from 55 days after the optical 
maximum through 251 days with the medium-energy pro- 
portional counter array (ME) and the low-energy imaging 
telescope (LE) onboard EXOSAT. They made six separate 
observations, 55, 62, 74, 83, 93, and 251 days after the op- 
tical maximum. The source became very weak in the ME 
from day 83 and only upper limits were derived. How- 
ever, there was still some residual X-ray emission detected 
even more than 250 days after maximum only with the 
LE band. The authors interpreted the early X-ray emis- 
sion in terms of emissions by the pre-outburst M-giant 's 
wind that had been shock-heated by the outburst ejecta 
(see also Bode & Kahn 1985; O'Brien & Kahn 1987; Itoh 
& Hachisu 1990; O'Brien, Bode, & Kahn 1992). On the 
other hand, the late X-ray emissions have been interpreted 
in terms of the continued hydrogen shell burning on top 
of the white dwarf. This interpretation is not consistent 
with our model calculation because the steady hydrogen 
shell-burning ended long before the day 251, i.e., at the 
day ~ 110. We will discuss this point again in the next 
subsection in relation to the distance to RS Oph. 

6.3. The Distance to RS Ophiuchi 

The distance to RS Oph is determined to be 0.57 kpc 

from the dereddened UV flux (Snijders 1987a) as shown 
in Figures 11 and 12. We assume the accretion disk exists 
during the outburst. We cannot reproduce the UV light 
curve if the disk does not exist or its size is smaller than 
a < 0.003, i.e., < 0.5i?Q, as already shown in Hachisu & 
Kato (2000a). The UV light curves depend hardly on the 
disk parameters of a, (3, or v when a > 0.01, mainly be- 
cause the UV light is coming from the innermost part of 
the accretion disk. The distance depends weakly on the 
irradiation efhciency (?7dk), for example, d = 0.70 kpc for 
Vdk = 1-0 or d = 0.50 kpc for r]uK = 0.25. The red gi- 
ant does not contribute to the UV because its irradiated 
temperature is as low as ~ 4000 K. Thus, the distance 
estimation is rather robust. 

We have calculated the absorption of Ay = 11.09 — 
51og(570/10) = 2.3 for d = 0.57 kpc from our light curve 
fitting. This absorption of Ay = 2.3 is consistent with 
the color excess of E{B-V) = 0.73 (Snijders 1987b) be- 
cause of Ay = 3.1 E{B — V). This distance of ~ 0.6 kpc is 
not consistent with other estimations of ti = 1.6 kpc from 
the hydrogen column density (Hjellming ct al. 1986), 
d = 1290 pc (Harrison et al. 1993), and d = 1 kpc 
(Sekiguchi et al. 1990), both of which are determined 
from the K-band luminosity. In the following, therefore, 
we discuss extensively other distance estimations support- 
ing our short distance of 0.6 kpc to RS Oph. 

The VLBI map 77 days after the optical maximum (Tay- 
lor et al. 1989) shows an elongated structure with a long 
axis of 0'.'22 and a short axis of 0'.'08. If we assume the 
distance of d = 1.6 kpc to RS Oph, this fastest aver- 
age expansion velocity should be 4,000 km s~^. This is 
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not consistent with the initial velocity of ejecta ~ 3, 000 — 
4, 000 km s~^ because the ejecta is strongly decelerated by 
the circumbinary matter that is formed by the red giant 
wind in quiescence. If the distance is d = 0.6 kpc, on the 
other hand, the average expansion velocity is much slower 
(~ 1,400 km s~^), which is consistent with the initial ve- 
locity of the ejecta. Then, the slowest average expansion 
velocity is estimated to be ~ 500 km s~^, which is also 
consistent with the optical line width (e.g., ~ 490 km s~^ 
110 days after optical maxiniiuri, Snijdcrs 1987a). 

The outer edge of the circumbinary matter (CBM) that 
is formed by the red giant wind is about iicBM ^ 1 x 10^^ 
cm, since an upper limit of 20 km to the velocity of 
the red giant wind has been given by Gorbatskii (1972) 
based on spectra by Dufay et al. (1964). Here, we assmne 
that the ejecta of the previous 1967 outburst swept up the 
previously existing CBM far away. The soft X-ray rapidly 
fell down from 60 70 days after the optical maximum so 
that the shock had broken out of the circumbinary matter 
on that day (Mason et al. 1987; Itoh & Hachisu 1990; 
Lloyd et al. 1993). If the outer edge of the circumbinary 
matter corresponds to the elongated radius of the VLBI 
map on the day 77 (O'.'ll of half axis), we obtain the dis- 
tance to RS Oph, d < i?cBM/(A6l) = 2 X lO^i cm = 0.65 
kpc. 

Very soft X-rays were observed 251 days after the optical 
maximum in the 1985 outburst of RS Oph (Mason et al. 
1987). Mason et al. estimated a blackbody temperature 
of 3.5 X 10^ K and a total energy flux of 

ix,obs = 1 X 10^7 ergss-\ (57) 

for a hydrogen column density of A^h = 3 x 10^^ cm~^, and 
concluded that these soft X-rays come from a WD pho- 
tosphere with steady hydrogen shell burning. However, 
steady hydrogen shell burning has stopped 109 days af- 
ter the optical maximum (HJD 2,446,202) in our model of 
X = 0.70, Z = 0.004, and Mwd = 1.377Mq (Figs. 11 and 
12), that is, this soft X-ray detection is indicating accretion 
luminosity rather than hydrogen shell burning luminosity. 
Using a blackbody temperature of Tbb ^ 3.5 x 10^ K and 
a WD radius of 0.003 Rq for 1.377 Mq, we estimate a 
total luminosity of 

47rR^j,aTiji ~ 120^0. (58) 

Thus, a rather short distance of 0.4 kpc to RS Oph 
can be derived by equating equation (57) and (58), i.e., 
Lx,ohs = Lx,wD- Thus, our short distance of 0.6 kpc is 
roughly consistent with the soft X-ray observation by Ma- 
son et al. (1987) about eight months after the optical 
maximum. This low luminosity corresponds to the accre- 
tion luminosity of 

M^c = ^-^^7;^^-^ ^ 0.2 X 10-^Me yr-^ (59) 

This accretion rate is six times lower than our estimated 
mass accretion rate of Mace = 1-2 x 1O~^M0 yr~^ averaged 
between 1967 and 1985. 

Dobrzycka et al. (1996) pointed out that the line 
fluxes in H I and He I before the 1985 outburst are 
about four times larger than the fluxes after the 1985 
outburst, thus indicating a decrease in the mass accre- 
tion rate after the outburst (see also the low luminosity 



during the post-outburst minimum ~100 — 400 days after 
maximum pointed by Evans et al. 1988). If our value 
of Mace = 1-2 x 1O~^M0 yr"'^ is adopted as the mass 
accretion rate before the 1985 outburst, these line fluxes 
before/after the outburst are roughly consistent with the 
ratio of the estimated mass accretion rates before and af- 
ter the outburst (1.2/0.2) if the line fluxes depend pro- 
portionally on the gravitational energy release. The day 
251 was just in the period of the post-outburst minimum 
{niv ~ 11.5) of the day ~100 — 400, after which the visual 
luminosity increased to my ^ 10.5, by about one magni- 
tude (e.g.. Fig. 1 of Evans et al. 1988). 

Bohigas ct al. (1989) estimated the H/3 intensity pho- 
toionized by X-ray photons from the central hot WD and 
compared it with the observation 201 days after maximum. 
Their theoretically estimated value of 

I{B.P) = 1.24 x 10"^"ergs cm^^s^^ 

is 3.5 times larger than the observed value of /(H/3) = 
3.55 X 10~^^ ergs cm~^ s~^ when we adopt i?,wD = 10^ cm 
and d= 1.6 kpc. However, if we use i?wD = 0.003 Rq and 
d = 0.6 kpc, the estimated value is /(H/j) = 3.8 x 10"-^^ 
ergs cm^^ s~^, being consistent with the observed inten- 
sity. 

X-rays were also observed at the quiescent phase (Orio 
1993), but the flux is too low to be compatible with the 
TNR model. One possible explanation is an absorption by 
the massive cool wind (e.g.. Shore et al. 1996) from the 
RG component as suggested by Anupama & Mikolajewska 
(1999). On the other hand, soft X-rays were observable 
just after the outburst, because the ejecta swept away the 
cool wind around the RS Oph system. The cool red giant 
wind fills again its Roche lobe and can cover the hot com- 
ponent in a time scale of Ai = a/u ~ 3OOi?0/lO km = 
300 days, which is roughly consistent with the recovery of 
the quiescent V luminosity as shown in Figure 1 of Evans 
et al. (1988). 

Hjellming et al. (1986) estimated the distance to 
RS Oph of d = 1.6 kpc from H I absorption line mea- 
surements, using the hydrogen column density of A^h = 
(2.4 ± 0.6) X 10^1 cm-2 and the relation of Nn/iT^ d) = 
1.59 xlO^" cm"2 j^-i ]jp^.-i together with Ts = 100 K. If a 
large part of this hydrogen column density stems not from 
the Galactic absorption but from the local one belonging 
to RS Oph, the distance is overestimated. For example, 
in the direction of the recurrent nova U Sco, Kahabka et 
al. (1999) reported a hydrogen column density of (3 — 
4)xl0^^ cm~^, which is much higher than the Galactic 
hydrogen column density of 1.4 x 10^^ cm~^. If RS Oph 
is in the same case as in U Sco, the distance of 1.6 kpc is 
just an upper limit and we possibly derive a much shorter 
distance. 

6.4. The Ejected Envelope Mass and Possibility of SN la 

Explosion 

The envelope mass at the optical maximum is estimated 
to be AM = 2.2 x 10"^Mq, which indicates a mass ac- 
cretion rate of 1.2 x 1O^^M0 yr^^ during the quiescent 
phase between the 1968 and 1985 outbursts. In our cal- 
culation, about 90% of the envelope mass (2.0 x 1O~^M0) 
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has been blown off in the optically thick wind, and the 
residual 10% (0.2 x lO-^Mg) has been left and added 
to the helium layer of the WD. The residual mass itself 
depends on both the hydrogen content X and the WD 
mass. It is roughly ranging from 10% {X = 0.70) to 
15%{X = 0.50) including ambiguities of the WD mass 
(MwD = 1.377 ± O.OIMq). Therefore, the net mass in- 
creasing rate of the WD is Mhc = 1-2 x 10"*^ Mq yr~^ 

Our calculated ejecta (wind) mass seems to be inconsis- 
tent with observation. Bohigas ct al. (1989) estimated 
the (fully ionized hydrogen) shell mass from their obser- 
vation 201 days after maximum. The obtained value is 
Msheii ^ 1-5 X lO-^{d/l kpcfMQ. Then, they obtained 
the ejecta mass form their equation (8), i.e., Mojccta ~ 
3{Vs/Vc)^Mshe\\, where Vs is the shell velocity after the 
shock has broken out of the circumstellar envelope and 
Ve is the velocity of the ejecta before it collides with the 
circumstellar matter. Assuming Vs = 200 km s"^ and 
Ve = 2000 km s~^, Bohigas et al. estimated Mejecta ~ 
4.5 X 10~^(d/l kpc)^M0. This ejecta mass is not consis- 
tent with our numerical model (AM^md = 2.0 x 1O~^M0) 
if we use d — 0.6 kpc. 

This inconsistency comes from the inappropriate pair 
of velocities because the estimated mass depends sensi- 
tively on the ratio Vg/Vg. If we adopt an appropriate pair 
of velocities, our wind model is consistent with the ob- 
servational shell mass. Starting from the set of d — 0.6 



kpc and M, 



ejecta 



AM, 



wind 



2.0 X 10"*^ Mq, we obtain 



Msheii ~ 5.4 X IQ-'^Mq and V^/V^ ~ 0.35. Adopting 
Ve ~ 4000 km s~^ from the observation (e.g., Rosino & 
lijima 1987; Snijdcrs 1987a,b) gives V; 1400 km s^i. 

This velocity of 1400 km s^^ is consistent with the ra- 
dio map obtained 77 days after maximum (Taylor et al. 
1989). The VLBI map shows an elongated structure of an- 
gular dimensions of ~ 0'.'22 (long axis) and ~ O'.'OS (short 
axis), which is consistent with the distance of 0.6 kpc and 
the fastest expansion velocity of 1400 km s~^ after shock, 
i.e., 0'.'092 {d/1 kpc)-i • (14/1000 km s"^) ~ 0'.'22. Thus, 
we have obtained a consistent set of physical quantities, 
Msheii ~ 5.4 X 10-6Mq, Mejccta ~ 2.o'x lO-^M©, and 
Vg = 1400 km s~^ for the distance of d = 0.6 kpc. 

Thus, we obtain the cool wind mass loss rate of 4 x 
1O~^M0 yr~^ (= Msheii/18 yr), indicating that one fourth 
of the cool wind has been captured by the WD. Obser- 
vationally, a very high mass loss rate of the cool wind, 
~ 10'"'' M3 yr^-'^, was suggested by Shore ct al. (1996) 
based on the UV data. However, Dobrzycka et al. (1996) 
argued against such a high mass loss rate of the M-giant 
cool wind. If the wind mass loss rate approaches such 
a high rate of IO^^Mq yr"-'^, a very large 10/im excess 
should be observed. But, the giant of RS Oph has only 
a modest 10 — 20/im excess. They estimated a few times 
10"'^ Mq yr^^ (e.g., Kcnyon, Fernandez-Castro, & Sten- 
cel 1988). This mass loss rate of the cool wind is very 
consistent with the present estimation. 

Finally, we examine whether or not RS Oph will explode 
as a Type la supernova in quite a near future. One of the 
conditions for SN la explosions is changed as follows: 

Mace > MHe-det ~ 4 X 10'^ Mq yr'S (61) 
for low metallicity stars (Nomoto 1982; Nomoto & Kondo 
1991) instead of equation (52). Our estimated mass accre- 
tion rate of Mace = 1-2 x 1O~^M0 yr~^ meets the above 



condition. Thus, we may conclude that RS Oph is an 
immediate progenitor of an SN la even if the donor is a 
metal-poor star. 

6.5. Absence of Secondary Maximum 

Am(mg the four recurrent novae with a red giant com- 
panion, only T CrB shows a secondary maximum of AV 
2 mag brighter than the quiescent level 120 — 190 days after 
the optical maximum. We discuss in more detail on the 
difference between T CrB and the others. The essential 
difference exists in the fact that the red giant component 
in T CrB fills its inner critical Roche lobe but the others 
do not (see the following sections for V745 Sco and V3890 
Sgr). 

In the T CrB system, the mass transfer continues from 
the red giant component to the white dwarf even in the 
outburst (hot) wind phase. Moreover, the strong irra- 
diation increases the mass transfer rate from the red gi- 
ant component; thus increasing the radius of the optically 
thick region of the accretion disk. The unstable warping 
grows shortly after the hot wind stops because the mass 
transfer rate satisfies the warping condition (49). 

On the other hand, the RS Oph system cannot main- 
tain the mass transfer because the cool red giant wind has 
once been evacuated by the outburst (hot) wind/ejecta 
from the white dwarf. During the day — 80 (hot wind 
phase), the growth of warping was suppressed by the hot 
wind because the hot wind (matter) momentum is much 
larger than the photon momentum as shown by Kato & 
Hachisu (1994). During the day 80 — 120 (steady burning 
phase), the warping instability may grow to some extent. 
However, the disk cannot contribute much to the V light, 
since the optically thick region of the accretion disk has 
become as small as or smaller than ^ l-R©- 

The mass transfer to the white dwarf resumes ^ 400 
days after the optical maximum, that is, after the hot wind 
stops (~ 80 days) and then the cool red giant wind fills 
again the binary system {At = a/v ~ 300 i?0/lO km s~^ 
— 300 days) as already discussed earlier. At that time, 
the luminosity of the white dwarf has once decreased as 
low as ^ 100 — 200 Lq and too low (two orders of mag- 
nitude lower) to excite the warping instability as easily 
understood from the condition for RS Oph (Southwell et 
al. 1997), i.e., 

1/2 

< 



3 X 10-8 Mq yr-i 



■disk 



i?WD 

0.003 Rq 



^bol 



2 X 10^8 ergs s-^ 



1/2 



MwD 

1.377 Mq 



-1/2 

(62) 



In the cool wind accretion phase, our estimated accretion 
rate of the WD in RS Oph is Mace 1.2 x IO-'^Mq yr"! 
(therefore, Lboi ~ 5 x 10^® ergs s^^ for the accretion lumi- 
nosity) , which also does not meet the above condition even 
if iidisk = 10 Rq, so that the radiation-induced instability 
does not occur in RS Oph. This is the reason why RS Oph 
does not show a second peak. 

The other two, V745 Sco and V3890 Sgr, are similar 
systems to RS Oph, so that a secondary maximum was 
not observed, as will be shown in the following sections. 

6.6. Summary of RS Oph Outburst 
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The main features of the RS Oph outbursts can be 
well understood with a thermonuclear runaway model on 
a very massive white dwarf close to the Chandrasekhar 
mass limit. We have summarized the following important 
points: 

(1) It has been suggested that the M-giant in RS Oph 
is metal-poor. Assuming an upper limit of white dwarf 
mass, i.e., Mwd = 1-377 Mq, we have obtain a lower limit 
of metallicity, i.e., Z > 0.004 from the light curve fitting of 
the early phase (~ 4 clays). In other words, the early 
visual light curve can be well reproduced only by a bloated 
white dwarf photosphere of Mwd = 1.377 ± 0.01 Mq 
for a low metallicity of Z = 0.004. On the other hand, 
Mwd = 1.35 ± 0.01 Mq is obtained for the solar metallic- 
ity ofZ = 0.02. 

(2) The ensuing phase (~ 5 — 80 days) visual light curve 
cannot be reproduced only by a white dwarf photosphere. 

It requires an irradiation effect of the accretion disk, the 
size of which is gradually decreasing from ~ 10 Rq to 
-li?0. 

(3) The UV light curve can also be well reproduced with 
the same model as that of the visual luminosity. The early 
phase UV flux mainly comes from the bloated white dwarf 
photosphere while the late phase UV flux comes from the 
irradiated disk surface because the main emitting region 
of the white dwarf photosphere has shift from UV to soft 
X-ray in the late stage. 

(4) The optically thick wind stops on the day ^80. The 
WD photosphere quickly shrinks. The sharp drop in the 
UV flux around the day 80 can be understood from this 
quick shrink of the WD photosphere. 

(5) The distance to RS Oph is estimated to be 0.6 kpc 
from the fitting with the dereddened UV flux. This dis- 
tance is consistent with the V light curve fitting. 

(6) The envelope mass of the bloated white dwarf is 
2.2 X 1O~^M0, suggesting the mass accretion rate of 
1.2 X lO-^M© yr-i. 

(7) The wind carries away about 90% of the envelope mass, 
so that the residual 10% can be accumulated to the helium 
layer of the white dwarf. As a result, the white dwarf grows 

at a rate of 1.2 x lO^^M© yr"^ 

(8) RS Oph will certainly explode as a Type la supernova 
in quite a near future if the white dwarf consists of carbon 
and oxygen. 

7. V745 SCORPII 

V745 Scorpii underwent the second recorded outburst 
in 1989 (the first recorded outburst in 1937), and has been 
recognized as a member of recurrent novae (e.g., Sekiguchi 
et al. 1990, and references therein). Its light curve is char- 
acterized by a very rapid rise (in about half a day) and 
decline {ts ^ 9 days, is the time taken to drop 3 mag 
from maximum). The general spectral evolution closely 
resembles that of RS Oph (e.g., Sekiguchi et al. 1990). 
Although its orbital period, the most important system 
parameter, is still not known, the secondary cool compo- 
nent has been classified as M6 — 8 III (from optical spectra, 
Sekiguchi et al. 1990; Duerbeck & Seitter 1990; WiUiams 
et al. 1991), M4 III (from the infrared spectra and colors, 
Harrison et al. 1993), or M4 — 6 III (from the absorption 
line indices, Anupama & Mikolajewska 1999). Therefore, 
V745 Sco is a twin system to RS Oph, a member of the 



recurrent novae with a red giant secondary. 

7.1. The Model of the 1989 Outburst 

We have calculated V light curves for five cases of the 
WD mass (Mwd = 1-3, 1.35, 1.36, 1.37, and 1.377Mq) 
by assuming the hydrogen content X = 0.70, 0.50, or 0.35 
with the solar metallicity {Z = 0.02). By fitting light 
curves in the first 25 days of the 1989 outburst, we have de- 
termined the WD mass of Mwd = 1.35±O.O1M0 as shown 
in Figure 13. Here, we have shifted the visual magnitude 
obtained with orange filter or 2415 films (Liller 1989a,b) 
by ~ 0.8 mag down (crosses, x) because they are sys- 
tematically about 0.8 mag brighter than the ^-magnitude 
obtained by Sekiguchi et al. (1990). The decline rate of 
the V light in the early phase depends very sensitively on 
the WD mass, but depends hardly on the hydrogen con- 
tent X, the companion (its mass and radius), or the size 
of the accretion disk, so that we can safely determine the 
WD mass from the early light curve (see also Kato 1999). 

Our fitting indicates an apparent distance modulus of 
(m - M)v = 17.0. Thus, the distance to V745 Sco is 
derived to be d = 5.0 kpc if we adopt the absorption of 
Av = 3.1 X E{B-V) = 3.5, where E{B-V) = 1.1 (Har- 
rison et al. 1993; Williams 1994). 

From our solutions, the envelope mass at the maximum 
V fight is estimated to be AM = 4.6 x IQ'^Mq, indicat- 
ing a mass accretion rate of Mace = 0.9 x lO^"^ Mq yr~^ 
between the 1937 and 1989 outbursts. About 95% of the 
envelope mass has been blown off by the wind and the 
residual 5% (2.3 x 10"'' Mq) has been left and added to 
the helium layer of the WD. Therefore, we obtain the net 
mass increasing rate of the WD as Mrb = 0.45 x 1O~^M0 

7.2. Irradiation Effect of Companion M-giant 

In the late phase {t ~ 25 — 200 days after the optical 
maximum) , contributions of the irradiated red giant (RG) 
and accretion disk (ACDK) to the V light may become 
dominant like as in T CrB and RS Oph. First, in this sub- 
section, we examine the irradiation effect of the M-giant, 
and then, in the next subsection, we also study the irradi- 
ation of the ACDK. 

In order to examine companion's irradiation effect, we 
need to know the separation, mass ratio, inclination an- 
gle, and the radius of the red giant, but none of them is 
known. The lack of both the orbital information and the 
visual luminosity observation in the late phase clouds the 
accurate determination of the outburst model because of 
ample ambiguity in the system parameters. 

The orbit is assumed to be circular. It is observationally 
suggested that orbital periods of S-type symbiotic stars are 
related to their spectral type by 

P ~ (1.31)^ x 190 days, (63) 
where the power of S means the sub-spectral type of M- 
giants (5 = 4 for M4 III, see, e.g., Miirset & Schmid 1999; 
Harries & Howarth 2000). From this relation, we assume 
P = 560 days for V745 Sco (M4 HI, e.g., Harrison et al. 
1993) and the ephemeris for inferior conjunction of the M- 
giant in front, (^min = HJD 2, 448, 700.0 + 560S. In what 
follows, we refer this ephemeris simply as (jiram = 700 day. 

It has also been suggested that M-giants in S-type sym- 
biotic stars lies well within the inner critical Roche lobe, 
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that is, 7 ~ 0.3 - 0.6 (e.g., Miirset & Schmid 1999). 
Therefore, we assume roughly the ratfo of 



In additfon, we also assume 50% efficiency of the irradia- 
tion (77RG = 0.5). For the companion M-giant, Sckiguchi 
et al. (1990) suggested my = 17 in quiescence from 
the V magnitude just after the outburst. Anupama & 
Mikolajewska (1999) obtained mssoo = 17.7 in March 
1998. We adopt here mv = 17.7 in quiescence. 

The light curves are calculated for six cases of the com- 
panion mass, i.e., Mrg = 0.5, 0.6, 0.7, 0.8, 1.0, and 
1.2Mq. Since we have obtained similar light curves for 
all of these six masses, wc show here only the results 
for Mrg = I.OMq. For this case, the separation is 
a = 38O.1-R0, the effective radii of the inner critical Roche 
lobes are i?J = 154. Oi?© for the WD component, and 
i?2 = 134.3i?0 for the RG component. If 7 = 0.5, then we 
have i?2 = -Rrg = 0.5i?2 ^ Q7Rq. Then, the nonirradi- 
ated photospheric temperature of the RG is estimated to 
be Tph,RG = 2700 K from the light curve fitting in quies- 
cence. 

We have finally adopted the best fit set of the inclina- 
tion angle i = 80° for the ephcmeris of (f>nun = 700 day 
as shown in Figure 13. Dash-dotted curves in Figure 14 
show V light curves for various inclination angles (i = 30°, 
60°, 80°) of the orbit for the ephemeris of ^min = 700 day. 
Their light curves depends hardly on the inclination angle 
because the non-irradiated hemisphere of the M-giant faces 
towards the Earth. On the other hand, V light is much 
brighter (dotted curves) if the ephemeris of (pmin = 500 
day is adopted because the irradiated hemisphere of the 
M-giant faces towards the Earth, which is not compatible 
with the observation. 

Figure 15 shows the dependency of the hydrogen con- 
tent, X. The early-phase V light curve depends hardly on 
the hydrogen content, while the late-phase V light curve is 
weakly identified by the hydrogen content of the WD en- 
velope. The optically thick, outburst wind stops 59, 69, 82 
days after maximum, and the steady hydrogen shell burn- 
ing ends 71, 90, 123 days after maximum, for X = 0.35, 
0.50, and 0.70, respectively. The V light curve indicates 
that the hydrogen content is somewhere between X = 0.5 
and X = 0.7 (dash-dotted curves). 

We have also calculated light curves for a much longer 
orbital period {P = 1000 days, for M6 III, e.g., Anupama 
& Mikolajewska 1999), and have obtained similar results 
to those for P = 560 days. In this case, the binary param- 
eters are a = 559.4i?Q, Rl = 226.7ii0, i?| = 197.6i?0, 
R,2 = 99i?0. We have estimated Trg = 2500 K for the 
best fit model in the late phase. 

7.3. Small Irradiation Effect of Accretion Disk 

In this subsection, we show that the ACDK should not 
contribute so much to the V light. Since the visual light 
curve is well reproduced by M-giant's irradiation with no 
accretion disk, one may imagine that the accretion disk 
does not contribute much. Here, we determine the upper 
limit of the ACDK size in relation to the reason why a 
secondary maximum does not occur in V745 Sco. 

We here again assume 50% efficiency of the ACDK irra- 
diation (r/DK = 0.5). The temperature of the unheated 



surface of the ACDK including the rim is assumed to 
be Tdisk = 1000 K. We have examined another case of 
^disk = K, but found no significant difi^erences in the V 
light curves. The power of v is assumed to be i/ = 9/8 
from the standard disk model. 

We have finally obtain allowable sets of the disk param- 
eters of a < 0.01 and (3 < 0.01 for the orbital inclination 
angle of i ~ 80° as shown in Figure 13. Here, the thin 
solid line denotes the case of a = 0.01 and (3 = 0.01, while 
the dash-dotted line corresponds to a = (no accretion 
disk). Both can reproduce the V light curve of the 1989 
outburst. 

We have examined six cases of the inclination angle of 
the orbit, i = 80°, 70°, 60°, 50°, 40°, and 30° and shown 
three of them in Figure 14. The V light curves are more 
luminous for a lower inclination angle of the orbit, since a 
larger area of the ACDK can be seen from the Earth. This 
means that lower inclinations of the orbit contradicts the 
existence of a large ACDK. To summarize, for any case, 
the ACDK should not contribute so much to the V light. 

Figure 15 shows the dependency of the hydrogen content 
for the case with an accretion disk. The V light curves sug- 
gest that the hydrogen content is not clearly identified by 
the observation but it may be somewhere between X = 0.5 
and 0.7 (solid curves). 

Thus, a very small size of the optically thick ACDK, i.e., 
^disk ^ 1 Rq, in the late phase indicates no secondary 
maximum of the light curve as discussed in RS Oph even 
if the warping instability sets in. 

7.4. The Distance, Envelope Mass of the White Dwarf, 
and Possibility of SN la Explosion 

The distance to V745 Sco is estimated to be d = 5.0 
kpc with the absorption of Ay = 3.5. This is roughly con- 
sistent with Harrison et al.'s (1993) estimation of d = 4.6 
kpc. Harrison et al. attributed M4 HI to the cool com- 
ponent from the infrared CO lines. Assuming the abso- 
lute magnitude of Mk = — 5.5 (M4 III), they obtain the 
distance modulus of (to — M)o = niK — Mk — Ak = 
8.2 - (-5.5) - 0.4 = 13.3, where they adopted tok = 8.2 
from Sekiguchi et al. (1990) and used a relation between 
the absorptions of Ak and Ay given by Rieke & Lebofsky 
(1985), i.e., Ak/Av = 0.112. 

Sekiguchi et al.'s (1990) K-band luminosity tuk = 8.21 
(87 days after maximum) is about 0.2 mag brighter than 
Harrison et al.'s (1993) K-band luminosity tuk — 8.41 
(two years after the 1989 outburst). This is consistent 
with our models as seen in Figure 13, because the cool 
component had been strongly irradiated by the WD photo- 
sphere always during Sekiguchi et al.'s observations. Using 
rriK = 8.41, we obtain a distance modulus of (to — M)q = 
tuk-Mk-Ak = 8.41 - (-5.5) - 0.4 = 13.5, i.e., d = 5.0 
kpc, being coincided with our distance estimation. 

At the end of this section, we discuss the possibility 
of SN la explosion. We have already obtained the average 
mass accretion rate of Mace = 0.9 x 10"'' M© yr~^ between 
the 1937 and 1989 outbursts from the envelope mass of 
AM = 4.6 X IQ-'^Mq. This meets condition (52) for SN 
la explosions (Nomoto & Kondo 1991), For all the three 
cases of X = 0.70, 0.50, and 0.35, about 95% of the enve- 
lope mass has been blown off in the outburst wind and the 
residual 5% (2.3 x 1O~'^M0) has been left and added to the 
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helium layer of the WD. Therefore, we obtain the net mass 
increasing rate of the WD of Mrb = 0.45 x 10~^Mq yr~^. 
Thus, we may conclude that V745 Sco is also an immedi- 
ate progenitor of Type la supernova if the hot component 
is a carbon-oxygen white dwarf. 

8. V3890 SAGITTARII 

V3890 Sagittarii underwent the second recorded out- 
burst in 1990 (since the 1962 first recorded outburst), and 
has been recognized as a member of the recurrent novae 
(Jones 1990a; Lillcr 1990). Its light curve is character- 
ized by a very rapid rise (within two days) and fast decline 
(ts ~ 17 days; see, e.g.. Fig. 1 in Gonzalez-Riestra 1992). 
The general spectral evolution closely resembles that of 
RS Oph (e.g., Sckiguchi 1990). The cool component has 
been classified from optical spectra as M4 III (Sekiguchi 
1990), M8 III (Williams et al. 1991), or M5 III (Anupama 
& Mikolajewska 1999). The infrared spectra and colors 
indicate a spectral type M5 III (Harrison et al. 1993). 
Thus, V3890 Sgr is a similar system to RS Oph and V745 
Sco, a member of the recurrent novae with an M-giant 
secondary. 

8.1. The Model of the 1990 Outburst 

We have calculated the total V light of the WD pho- 
tosphere and the non-irradiated RG photosphere for the 
WD mass of Mwd = 1-2, 1.3, 1.34, 1.35, 1.36, 1.37, and 
1.377M0 by assuming the hydrogen content X = 0.70 and 
the solar metallicity of Z = 0.02. Four of them (1.2, 1.3, 
1.35, and 1.37 Af©) arc plotted in Figure 16 (solid curves). 
Here, we assume that the optical maximum was reached 
on HJD 2,2448,009. This is because the irradiation of the 
M-giant cannot contribute to the early phase V light. 

The observed V light curve in the first 20 days can be 
roughly reproduced by the model of 1.3M0 WD, but it 
deviates greatly in the ensuing phase (20 — 40 days after 
maximum). The envelope mass of the 1.3Mq WD at the 
optical maximum is too large to be compatible with the 
TNR model. The envelope mass of AM ~ 5 x 10~^Mq 
suggests Mace ~ 2 X 1O~''M0 yr"-'^, which is just below the 
lower limit of steady hydrogen shell burning for X = 0.7, 
i.e., Mstd in equation (7) and also in Figure 2. 

Only a bloated WD photosphere cannot reproduce the V 
light curve even in the early phase of outburst. Moreover, 
our simple picture of the WD photosphere never shows a 
sharp spike on the light curve as seen in Figure 16. These 
features strongly suggest that the irradiation of the ac- 
cretion disk is important even in the early phase: if the 
optically thick region of the accretion disk is as large as 
or larger than ~ 10 Rq , the decay of the V light becomes 
much slower as seen in RS Ohp (see, e.g.. Fig. 11). There- 
fore, we seek a reasonable set of system parameters based 
on the assumption that the V3890 Sgr system consists of 
the WD photosphere, the irradiated M-giant, and the ir- 
radiated accretion disk. 

8.2. Large Size of Accretion Disk in Early Phase 

In order to obtain a consistent set of system parame- 
ters, we have to adopt the WD mass more massive than 
M-^jffu ^ 1.35M0, and a size of the irradiated ACDK larger 
than iidisk ^ 10 Rq at the initial phase of outburst. The 



condition of Mwd ^ 1.35M0 is required from the consis- 
tencies with the TNR model, i.e., the mass accretion rate 



is smaller than Mace 5, 1 x 10"'^M( 



yr" 



and the rapid 



decline of the V light curve 20 40 days after maximum as 
shown in Figure 16. The V hght curve of a 1.34M0 WD 
model lies above the observational points on the day 26 
and the upper limits on the day 33 and 37, being not con- 
sistent with the observation. On the other hand, V light 
curve of a 1.35M0 or 1.37M0 WD is less bright during the 
day 3 — 40. The lack of the luminosity during these days 
can be supplied by the irradiation of a large size accretion 
disk. 

Figure 17 shows V light curves of fixed-size disk models 
for the WD masses of 1.35 Mq (solid curves) and 1.37 Mq 
(dashed curves). The disk size parameters are a = 0.7, 0.1, 
and 0.01. It is clear that the fixed-size disk model cannot 
reproduce the V light curve in the later stage. This is 
simply because the large disk contributes too much in the 
later phase. 

Here, we summarize our parameters: the orbital pe- 
riod of P = 730 days from equation (63) and the spectral 
type of M5 III (Anupama & Mikolajewska 1999) for the 
cool component, the ephemcris for the M-giant in front is 
?^min = HJD 2,448,050.0 -^730 xE,j = 0.5 from equation 
(64), 50% efficiency of the irradiation (t/rg = 0.5), the 
M-giant mass of Mrg = 1-0 Mq. In this case, the separa- 
tion is a = 453.5i?0, the effective radii of the inner critical 
Roche lobes are Rl = 183.8i?0 for the WD component 
and i?2 = 16O.2i?0 for the RG component. Then, we have 
Rrg = 0.5i?2 ~ 8Oi?0. The nonirradiated photospheric 
temperature of the RG is estimated to be Tph.Rc = 2600 
K from fitting with the quiescent V magnitude of my ^ 16 
(Anupama & Mikolajewska 1999). For the accretion disk, 
we asstune 50% efficiency of the irradiation {r]uK = 0.5), 
the unheated surface temperature of Idisk = 1000 K. 

8.3. Evaporating Accretion Disk 

The optically thick strong wind blows during the out- 
burst so that the disk surface is certainly blown off in 
the wind due to drag and ablation. Moreover, the mass 
transfer from the M-giant is certainly inhibited during the 
strong wind phase because the cool red giant wind has 
been evacuated by the hot wind from the WD as discussed 
in RS Oph. In such a situation, it is very likely that the 
edge of the disk is gradually dissipating with time. We 
here simulate such an effect by shrinking the disk size ex- 
ponentially, i.e., 

/q, \ (t-*o)/30 

a = ao[ — ] , (65) 

where we adopt = 0.7 and ai = 0.01, io is the time 
of the optical maximum, and the time t is in units of day. 
We set the minimum size of the ACDK at a = 0.001 in 
our calculations. 

Thus, we have obtained the best fit model for the WD 
mass of A'/wD = 1-35 ± 0.01 Mq as shown in Figures 17 
and 18. The envelope mass at the maximum V light is 
estimated to be AM = 3.1 x 1O~^M0, indicating a mass 
transfer rate of Mace = 1-1 x 1O~^M0 yr"-'^ between the 
1962 and 1990 outbursts. The efficiency of the ACDK irra- 
diation has been examined for other two cases of ?7dk = 1-0 
(dot-dashed line) and 0.25 (dotted line) only for X = 0.70 
and Mwd = 1-35 Mq as shown in Figure 18. 
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We have examined the dependency of the V hght curve 
on the orbital inchnation angle. The luminosity of the 
ACDK depends largely on the the orbital inclination an- 
gle mainly because a larger area of the ACDK can be seen 
from the Earth for a smaller inclination angle. Therefore, 
the inclination angle should be as low as « < 30° to re- 
produce the V light curve during t 5 — 20 days after 
maximum. 

We have also calculated V-light curves for other four 
cases of the companion mass, i.e.. A/rg =0.6, 0.7, 0.8, 
and 1.2 Mq, but obtained similar light curves to Mrg = 
1.0 Mq. We have checked the other case of Tdisk = 
K, but found no significant differences in the light curves. 
The power of v is fixed to be = 9/8, but we cannot find 
any significant differences in the light curves even if we 
adopt p = 2. Finally, we have study the effect of hydro- 
gen content. Either the early- or late- phase V light curve 
depends hardly on the hydrogen content X as shown in 
Figure 18. The optically thick, outburst wind stops 56, 
68, and 80 days after maximum, and the steady hydrogen 
shell burning ends 68, 89, and 122 days after maximum, 
for X = 0.35, 0.50, and 0.70, respectively. 

Anupama & Sethi (1994) obtained optical spectra 18 
20 days after maximum. They pointed out that these 
line features closely resembles the spectrum of RS Oph 
60 days after maximum. Anupama & Sethi derived the 
Zanstra temperature T« = 300, 000 K of the ionizing cen- 
tral source from the helium and hydrogen line (H/? and 
He II 4686 A) ratios and its radius i?* = O.O6ii0 from H/3 
line flux. For RS Oph, our model gives us Tph = 300, 000 
K and i?ph = 0.08 Rq about 70 days after maximum (see 
the previous section), being roughly consistent with the 
Zanstra temperature derived by Anupama & Sethi. For 
V3890 Sgr, however, our model suggests Tph = 83,000 
K and Rph = 1.0 Rq on the day 18, being not consistent. 
This may indicate the importance of the contribution from 
the shocked layer photons. 

Mukai et al. (1990) obtained optical spectrum on May 
14 (17 days after maximum). Comparing them with the 
earlier spectra observed by Wagner, Bertram, & Starrfield 
(1990) and Buckley et al. (1990), Mukai et al. pointed out 
that V3890 Sgr had progressed to a higher ionization state 
on that day. In our modeling of the light curve, the accre- 
tion disk plays a key role in these days. It is very likely 
that the conspicuous change in the spectrum ~ 20 days 
after maximum is caused by the evaporation of the accre- 
tion disk. As seen in Figure 17, the disk size decreases to 
< 1 Rq on the day ^ 20. Then, the mean temperature of 
the disk surface increases to > 50, 000 K. The same thing 
happens in RS Oph around the day ^^60 70. 

We thus suggest a very small size of the optically thick 
ACDK or no disk, that is, i?disk ^ 0.1 Rq, in the late 
phase (after HJD 2,448,036, indicated by "disk vanishes" 
in Figure 17). Therefore, no secondary maximum of the 
light curve is expected as well as in RS Oph and V745 Sco. 

8.4. The Distance, Envelope Mass of the White Dwarf, 
and Possibility of SN la Explosion 

Fitting indicates an apparent distance modulus of (m — 
M)v = 14.7. Thus, the distance to V3890 Sgr is de- 
rived to be d = 4.2 kpc if we adopt the absorption of 
Av = 3.1 X E{B-V) = 1.6, where E{B-V) = 0.5 



(Harrison ct al. 1993), because the distance modulus 
is (m - M)o = (m - M)v - Ay = 14.7 - 1.6 = 13.1. 
This is roughly consistent with Harrison et al.'s (1993) 
estimation of d = 5.2 kpc. Harrison et al. determined 
the infrared spectral type of the cool companion as M5 
HI. Assuming Mk = -5.5 (M5 III), they obtained the 
distance modulus of (m — M)o = niK — Mk — Ak = 
8.24 - (-5.5) - 0.17 = 13.6 by using Rieke & Lebof- 
sky (1985) relation of Ak/Av = 0.112. It should be 
noted here that there arc a large divergence in the ob- 
served color excesses: Gonzalez-Riestra (1992) deter- 
mined E{B-V) = 1.1 from UV lines; Williams (1994) 
obtained E{B — V) = 0.76 from H I and He II line ratios; 
Anupama & Sethi (1994) obtained E{B-V) = 1.1 from 
the He I line ratios, and E{B-V) = 1.2 from the Ha/H/3 
ratio. 

Finally, we examine the possibility of Type la supernova 
explosion. The envelope mass at the maximum V light is 
estimated to be AM = 3.1 x 1O~^M0 from our solutions, 
indicating a mass transfer rate of M^cc = 1-1 x 10~^Mq 
yr"-'^ between the 1962 and 1990 outbursts. Moreover, 
about 90% of the envelope mass has been blown off in 
the outburst wind and the residual 10% (3 x 10"'' M©) has 
been left and added to the helium layer of the WD. There- 
fore, we obtain the net mass increasing rate of Mhc = 
1.1 x IQ-^Mq yr-i. This meets condition (52) of SN la 
explosion (e.g., Nomoto & Kondo 1991), if the WD con- 
sists of carbon and oxygen. Thus, we may conclude that 
V3890 Sgr is an immediate progenitor of Type la super- 
nova if the hot component is a carbon-oxygen white dwarf. 
The white dwarf can grow to 1.378M0 even if mild helium 
shell flashes repeatedly occur after the helium layer grows 
to the critical mass at which helium ignites (see, e.g., Kato 
& Hachisu 1999). 

9. CONCLUSIONS 

We have developed a numerical method for calculating 
light curves of recurrent nova outbursts. Our basic model 
consists of a white dwarf with an accretion disk and a 
red giant companion. The calculated light curves includes 
contributions from the white dwarf photosphere, the irra- 
diated accretion disk, and the irradiated red giant. We 
summarize our main results on the light curve analysis of 
the four recurrent novae with a red giant companion, T 
CrB, RS Oph, V745 Sco, and V3890 Sgr. 

9.1. Relevance to Type la Supernova Progenitors 

1. The mass of white dwarf (WD) has been estimated, 
from the decline of early-phase light curve, to be Mwd = 
1.37±0.01 Mq for T CrB, 1.35±0.01 Mq with solar metal- 
hcity {Z = 0.02) or 1.377 ± 0.01 Mq with low metalhcity 
{Z = 0.004) for RS Oph, 1.35 ± 0.01 Mq for V745 Sco, 
and 1.35 ± 0.01 Mq for V3890 Sgr. 

2. The apparent distance moduli are obtained, with the 
light curve fitting, to be (m - M)v = 10.2 for T CrB, 
(m - M)v = 11.2 for RS Oph, (m - M)v = 17.0 for 
V745 Sco, and (m - M)y = 14.7 for V3890 Sgr. There- 
fore, the distances are estimated to be d = 0.94 kpc with 
Av = 0.35 for T CrB, d = 0.57 kpc with Ay = 2.3 for 
RS Oph, d = 5.0 kpc with Ay = 3.5 for V745 Sco, and 
d = 4.2 kpc with = 1.6 for V3890 Sgr. 
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3. Each envelope mass at the optical maximum is also cal- 
culated to be AM - 3 X IQ-'^M© (T CrB), 2 x IO'^Mq 
(RS Oph), 5 X IO-^Mq (V745 See), 3 x lO-^M© (V3890 
Sgr), which indicates an average mass accretion rate of 
Mace ~ 0.4 X lO-'^M© yr-1 (80 yrs, T CrB), 1.2 x IO-'^Mq 
yr-i (18 yrs, RS Oph), 0.9 x IO'^Mq yr"! (52 yrs, V745 
Sco), 1.1 X lO-'^Mg yr-i (28 yrs, V3890 Sgr) during qui- 
escent phase. 

4. Although a large part of the envelope mass is blown 
in the wind, eac;h white dwarf can retain a substantial 
part of the envelope mass after hydrogen burning. Thus, 
we have obtained net mass-increasing rates of the WDs, 
A/hc ^ 0.10 X IQ-HIq yr-i (T CrB), 0.12 x lO-^Af© yr"! 
(RS Oph), 0.05 X 10-^Mq yr"! (V745 Sco), 0.11 x IO'^Mq 
yr-i (V3890 Sgr). 

5. These results strongly indicate that the white dwarfs 
in recurrent novae have now grown up to near the Chan- 
drasekhar mass limit and will soon explode as a Type la 
supernova if the white dwarf consists of carbon and oxy- 
gen. 

9.2. Appearance/Absence of Secondary Maximum 

6. Among the four recurrent novae with a red giant com- 
panion, only T CrB shows a secondary maximum 120 — 190 
days after the optical maximum. This secondary max- 
imum can be well reproduced if an irradiated accretion 
disk of i?disk ~ 6 Rq around the white dwarf tilts due to 
a radiation-induced warping instability. If we neglect the 
accretion disk, we cannot fully reproduce the secondary 
peak by an irradiated M-giant model. 

7. The difference between T CrB and the others is in the 
companion: it fills its inner critical Roche lobe only in T 
CrB but in the others does not. Therefore, in T CrB the 
mass transfer is mainly driven by the Roche lobe overflow 
while in the other systems matter is fed by cool red giant 
winds. In T CrB, the mass transfer continues even in the 
outburst wind phase, which maintains an optically thick 



accretion disk large enough to show a secondary maximum 
after a tilting instability grows up. 

8. The other three systems, RS Oph, V745 Sco, and V3890 
Sgr, cannot maintain the mass transfer during the outburst 
because the cool red giant wind has once been evacuated 
by the ejecta of the outburst. The light curves of RS Oph 
and V3890 Sgr require a relatively large size accretion disk 
of i?disk ~ 10 Rq in the mid-phase but a small size of 
-Rdisk ^ 1 in the late-phase. The late-phase UV data 
of RS Oph can be well reproduced with a small accretion 
disk of i?disk ~ 0.5i?o. No or a very small accretion disk is 
also derived in V745 Sco. These accretion disks are quickly 
dissipating to become as small as i?disk ^ 1 Rq, at least, 
in the late phase of the outbursts, because no mass is fed. 
Such a small disk cannot contribute much to the V light. 
Even if the tilting instability begins to grow after the wind 
stops, it is therefore very unlikely that secondary maxima 
like in T CrB are observed in these three systems. 

9. The post-outburst minima of RS Oph outbursts, which 
is 1 mag (in V) below the quiescent level and continues 
for ^ 300 days from 100 to 400 days after the optical 
maximum, is understood from the evacuation of the cool 
red giant wind by the ejecta of nova. After the hot WD 
wind stops 80 days after the optical maximum, the cool 
red giant wind expands to recover again in 300 days. 
Then, ~ 400 days after maximum, mass accretion from 
the red giant wind to the white dwarf resumes. Supersoft 
X-rays from the white dwarf, whic;li were observed about 
250 days after maximum, are very consistent with the end 
of hot WD wind and the evacuation of the cool red giant 
wind. 
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Table 2 

Physical properties of Recurrent novae outbursts 
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^taken from Hachisu & Kato (2000a) of solar metallicity model for RS Oph 
^the present low metallicity {Z = 0.004) model for RS Oph 

*^taken from Hachisu et al. (2000a) 
'^taken from Hachisu & Kato (2000b) 
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Optically Thick Winds 




Fig. 1. — Optically thick winds blow from mass-accreting white dwarfs {left figure) when the mass transfer rate from a lobe-filling companion 
{right figure) exceeds a critical rate, i.e., Mace > Afcr- Here, we assume that the white dwarf accretes mass from equatorial region and, at the 
same time, blows winds from polar regions as illustrated in the figure. 
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Response of WD to accretion (X=0.70,Y=0.28,Z=0.02) 



-4 - 







■8 - 
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Optically Thick Wind 
(Accretion Wind) 



Steady H— burning 



(dM/dt) 

4: 




log(AM/MQ)= -3.5 



-12 



0.8 



1.2 



1.4 



Fig. 2. — Response of white dwarfs to mass accretion is illustrated in the white dwarf mass and the mass accretion rate plane, i.e., in 
MwD-Afacc plane. Strong optically thick winds blow above the line of Afacc > Mcr- The wind mass loss rate is M„ind ~ Afacc — Mct- 
Steady hydrogen shell burning with no optically thick winds occur between Mgtd ^ Mace < Mcr- There is no steady-state burning below 
Afacc < Afgtjj. Instead, intermittent shell flashes occur. The envelope mass at which hydrogen shell flash occurs is also shown (taken from 
Fig. 9 of Nomoto 1982). 
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final 
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Fig. 3. — Regions producing SNc la are plotted in the logP — (orbital period — donor mass) plane for the WD+MS system (left in 
the figure) and the WD+RG system (right in the figure). Here we assume an initial white dwarf mass of I.OMq. The initial system inside 
the region encircled by a thin solid line (labelled by "initial") is increasing its white dwarf mass up to the critical mass (Mia = 1.378Mq) 
for Type la supernova explosion, the regions of which are encircled by a thick solid line (labelled by "final"). Currently known positions of 
each recurrent nova are indicated by a star mark (*) for U Sco (e.g., Hachisu et al. 2000a), a filled triangle for T CrB (e.g., Belczyhski & 
Mikolajewska 1998), an open rectangle for T Pyx (e.g., Patterson et al. 1998), but by arrows for the other three recurrent novae, V394 CrA, 
CI Aql, and RS Oph, because the mass of the companion is not yet available explicitly. Two subclasses of the recurrent novae, U Sco subclass 
and RS Oph (or T CrB) subclass correspond to the helium-rich supersoft X-ray source channel (Hachisu et al. 1999b) and the symbiotic 
channel (Hachisu et al. 1999a) of Type la supernovae, respectively. 
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Fig. 4. — Wind mass loss rate (dashed lines) and total mass decreasing rate of the envelope (solid lines labeled by AM), i.e., wind mass loss 
rate plus nuclear burning rate, M^inj -f Mnuc, are plotted against the envelope mass for WDs with masses of 0.8, 1.0, 1.2, 1.3 and 1.377M0. 
There exist only wind solutions above the breaks on the solid lines while only static (no wind) solutions exist below the breaks. Thus, an 
optically thick wind stops at this break. The photospheric radius, Rph, the photospheric temperature, Tph, and the photospheric velocity, 
iiph, are also plotted against the total decreasing rate of the envelope mass in the figure. We have shown four different cases of hydrogen 
content of the envelope material, i.e., (a) X = 0.35, (b) 0.15, (c) 0.10, and (d) 0.05, from top to bottom panel, for the solar metallicity of 
Z = 0.02. Each sequence is plotted from the wind solution with log Tph = 4.0 to the static solution with Mnuc = l^std^ i-^-, the lower limit 
for steady hydrogen burning. 
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Fig. 5. — Same as Fig. 4, but for low metallicities oi Z = 0.01 (upper two panels) and Z = 0.004 (lower two panels), i.e., (a) X = 0.35, 
Z = 0.01, (b) X = 0.10, Z = 0.01, (c) X = 0.35, Z = 0.004, and (d) X = 0.10, Z = 0.004. 
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Fig. 6. — Model configuration near the second peak of the recurrent nova T CrB. The cool component {right in the figure) is a red giant 
filling up its inner critical Roche lobe. The hemisphere is heated up by the hot component (1.37 Mq white dwarf, left in the figure). The 
photospheric radius of the hot component near the second peak is as small as ~ 0.003 Rq, about ~ 0.00003 times the size of the cool 
component, but it is exaggerated in this figure to easily see it. We assume a warping accretion disk around the hot component, which is 
precessing at a period of about 140 days. The surface of the accretion disk is also heated up by the hot component. The radius of the 
accretion disk near the second peak is as small as ~ 6 -R0, i.e., about ~ 0.07 times the inner critical Roche lobe size. The accretion disk is 
also exaggerated in this figure to easily see it. 
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B 



Fic;. 7. - Each patch is divided into four triangle parts as shown in the figure, because each patch is almost a rectangular, but strictly 
speaking, not a rectangular for the Roche geometry of the companion star and for the warping disk surface. So that the normal unit vector 
to each patch, n, is approximately calculated from the four triangle parts by equation (22). 
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Fig. 8. — When patch k blocks the light from patch j to patch i, we do not sum up the contribution from j to i as shown in equations (29), 
(30), and (31). 
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Fig. 9. — Model light curves are plotted against time (HJD 2,430,000+) together with the observational points {open circles: Pettit 
1946a,b,c,d). Dotted: the total V light of the white dwarf (WD) photosphere and the red giant (RG) photosphere without irradiation. Large 
arrows attached by "min" indicate epochs at the spectroscopic conjunction with the M-giant in front. An ellipsoidal variation can be seen in 
the late phase. Dash-dotted: the total V light of the WD photosphere and the RG photosphere irradiated by the WD. The observed second 
peak cannot be reproduced only by the irradiation of the M-giant. Solid: the total V light of the WD photosphere, the RG photosphere with 
irradiation, and the accretion disk (ACDK) surface heated-up by the hot component. The parameters specifying the light curves are shown 
in the figure. 
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RS Oph 1985 Outburst (m-M)y=11.2 
X=0.70, M^r)=1.377 Mq WD photosphere only 




HJD 2,446,000+ (day) 

Fig. 10. — Model V light curves for various metallicitics are plotted against time (HJD 2,446,000+) together with the observational points 
of RS Oph outbursts. Filled circles indicate observational points with the previous outbursts superposed (taken from Rosino 1987, see also 
Haehisu & Kato 2000a). Open circles correspond to the observational points of the 1985 outburst, which are taken from lAU Circulars (Scovil 
1985; Lowdcr 1985; Alcock 1985; Hurst 1985a,b; Danskin 1985; Ducoty 1985a,b,c; McNaught 1985a,b,c,d; Lubbock 1985a,b; King 1985; 
Weier 1985; Riggs 1985; Cook 1985a,b,c; Mcdway 1985; Bortle 1985a,b,c; EUerbe 1985; Verdenct 1985; Sventek 1985; Albinson et al. 
1985), i.e., No.4031, 4036, 4048, 4049, 4066, 4074, and 4091. The WD mass is assumed to be the extreme case of AfwD = 1-377 Mq just 
before the SN la explosion. The metallicity is attached to each line (for Z = 0.001, Z = 0.004, Z = 0.01, Z = 0.02). Each line denotes the 
V light of the WD photosphere only. The hydrogen content of the WD envelope is assumed to be X = 0.70 for all models. The apparent 
distance modulus is (m — M)y = 11.2. The lowest metallicity case of Z = 0.001 cannot reproduce the rapid decline of the V light curve of 
RS Oph outburst. Since the hydrogen content, X, hardly affect the decline rate of the V light curve (e.g., Kato 1999), we may conclude that 
Z > 0.004 in RS Oph outbursts. 
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Fig. 11. — Same as Fig. 10 but for the entire outburst phase for MwD = 1.377M0, X = 0.7, and Z = 0.004. Dashed line: Sum of the 
contributions from the white dwarf (WD) photosphere and the non-irradiated red giant (RG) photosphere with the surface temperature of 
Tph = 3400 K. The RG lies well within the inner critical Roche lobe, i.e., its radius is 0.25 times the Roche lobe size (7 = 0.25). Dashed- 
doted: Sum of the WD photosphere and the irradiated RG photosphere. Solid: the total V light of the WD photosphere, the irradiated 
RG photosphere and the accretion disk (ACDK) surface. Dotted: ultra-violet (UV) radiation (911-3250 A) from the WD jjliotosphere, the 
irradiated RG photosphere, and the irradiated ACDK surface. The UV data (large open circles) arc taken from Siiijdcrs (1987a). Here both 
the disk shape and size are assumed to be constant during the outburst, i.e., a = 0.01 (-Rdisk ~ 14 Rq) and ,8 = 0.05 for the upper solid 
(and upper dotted) line, and a = 0.008 (i?disk ~ 1 ^0) ^^d f3 = 0.05 for the lower solid (and lower dotted) line, respectively. Two epochs 
with "wind stops" and "end of steady H burning" are indicated by arrows. The optically thick wind blows during the period from the first 
phase of the outburst (HJD 2,446,092) to 79 days after maximum (HJD 2,446,170). The steady hydrogen shell burning ends at 112 days after 
maximum (HJD 2,446,203). 
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Fig. 12. Same as Fig. 11 but for a case of gradually evaporating accretion disk. Here, we assume that the radius of the accretion disk is 
gradually shrinking from ~ 14i?Q {a = 0.1) to ~ IRq {a = 0.008) for X = 0.7, to ~ 2Rq (a = 0.012) for X = 0.5, to ~ 3Rq {a = 0.02) for 
X = 0.35, during the wind phase. Its time variation is exponential as given by equation (56). 
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Fig. 13. — Model light curves axe plotted against time (HJD 2,447,000+) together with the observational points of the V745 Sco 1989 
outburst. Here, the inclination angle of i = 80° is assumed. The other system parameters are shown in the figure. Filled circles indicate 
the observational points with V filter (taken from Sekiguchi et al. 1990). Crosses (taken from Liller 1989a,b) are shifted by 0.8 mag down 
because they are systematically 0.8 mag brighter than the V magnitude. Small allows (downwards) indicate an upper Umit before the 1989 
outburst (taken from Honda 1989). Open circles are visual magnitudes taken from Shore et al. (1989) and Overbeek (1989). Dashed line 
denotes the V light curve of the WD photosphere and the nonirradiated RG only. Dash-dotted line corresponds to the V light from the WD 
photosphere and the irradiated RG. Solid line represents the total V light of the WD photosphere, the irradiated RG, and the irradiated 
accretion disk. The photospheric radius of the WD is indicated above the V light curve. Two epochs are indicated by arrows with "wind 
stops" and "end of steady H-burning." The optically thick wind stops at HJD 2,447,819. The steady hydrogen shell burning ends at HJD 
2,447,860. Fitting with the observations indicates an apparent distance modulus of (m — M)v = 17.0. 
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Fig. 14. — Same as Fig. 13 but for different inclination angles (i = 30, 60, 80° from top to bottom). Dot-dashed lines denote the model light 
curves for the cases of no accretion disk while solid lines denote the cases with an accretion disk. The model light curves with an accretion 
disk depend largely on the inclination angle, while we do not find any significant differences for the cases of no accretion disk. Dotted lines 
indicate a different ephemeris of (pmin = 500 day, which seems not to be compatible with the observation for all three inclination angles. 
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Fig. 15. — Same as Fig. 13 but for three different hydrogen content oi X = 0.35, 0.50, and 0.70. We do not determine the hydrogen content 
from these observational points, but may conclude that it is somewhere between X = 0.5 and 0.7. 
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Fig. 16. - Model light curves are plotted against time (HJD 2,448,000+) together with the observational points of the V3890 Sgr 1990 
outburst. Here, the inclination angle of i = 30° is assumed. The other system parameters arc shown in the figure. Open circles indicate the 
observational points taken from lAU Circulars (Jones 1990a, b,c; Pcarce 1990a, b,c,d,c,f; Pereira 1990a, b; Boattini 1990: Lillcr 1990; Schmeer 
1990a, b,c) and Figure 1 of Gonzalez- Ricstra (1992), while small downward arrows show an upper limit. Solid lines denote the total V light 
of the WD photosphere and the nonirradiatcd M-giant secondary. Dot-dashed lines depict the total V light of the WD photosphere arid the 
irradiated M-giant secondary with two different </>mini i-c, ipmin = 50 day (lower) and (pmin = 250 day (upper) only for AfwD = 1-35 M3. 
Fitting with the observations suggests an apparent distance modulus of {m — M)y = 14.7. Here, we assume that the WD photosphere reached 
its maximum expansion on HJD 2,448,009. 
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Fig. 17. — Same as Fig. 16 but for two white dwarf masses, i.e., 1.35 Mq (solid lines) and 1.37 Mq (dashed lines). Both solid and dashed 
lines represent the total V liglit of the WD photosphere, the irradiated M-giant, and the irradiated accretion disk, the size of which is fixed 
during the outburst as defined by equations (11) and (12). The optically thick wind blows during the outburst, from the beginning of the 
outburst (HJD 2,448,008) to 80 days after maximum (HJD 2,448,089), and then the steady hydrogen shell-burning ends at 122 days after 
maximum (HJD 2,448,131) for the 1.35 Mq WD and the hydrogen content of X = 0.7, while 57 days and 78 days, respectively, for the 1.37 
Mq WD and the hydrogen content oi X = 0.7. The slow decline just after the optical maximum can be reproduced by the relatively large 
disk size around the 1.35 Mq WD but not by the 1.37 M© WD. 
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Fig. 18. — Same as Fig. 17 but for the case of evaporating accretion disk with three different hydrogen content of X = 0.35, 0.50, and 0.70 
{solid lines). We do not pose any constraints on the hydrogen content from these observational points, because of a lack of data points in the 
later phase. For comparison, two other cases with different irradiation efficiencies, ?7dk = 1-0 (dot-dashed) and 0.25 (dotted), are added only 
for X = 0.70. 



